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ABSTRACT 


An analysis has been carried out of a data base of polarized long wave infrared images 
of the instrumented Research Vessel POINT SUR recorded over a period of two days during 
the EOPACE measurements series in San Diego Bay in 1996. The measurements were made 
from a land site on Point Loma with an AGA780 sensor with internally mounted polarization 
filters. The objectives of the analysis were to determine a possible influence of target aspect 
- angle on the polarization signature, to compare polarization contrast improvement in San 
Diego Bay with previous measurements in the North Atlantic, and to validate by measurement 
the estimation of unpolarized signature from vertical and honzontal components. 5508 images 
representing 70 cases with vertical, horizontal and unpolarized sequences were analyzed. 
Using a horizontal polarizer, target to background contrast improvement was found with a 
mean of 1.08 (8%) compared with the 15% found in previous measurements. Estimated 
unpolarized signatures from vertical and horizontal components agreed with unpolarized 
measurements with a slope coefficient of .85 to .99. Target signature for major ship facets and 
for total ship showed no discernable degree of polarization. A total of 37 IDL programs 


developed for this analysis can be assembled as a package for future data processing. 
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I. INTRODUCTION 


Infrared systems have been shown to be among the most important sensors for the 
future sea warfare scenario. Intensive scientific and technological efforts have been made 
to improve the performance of those systems. The characteristic operation in a very 
cluttered environment enforces a requirement to improve the target-to-background 
contrast to improve discrimination. Analysis of polarization effects in the infrared target 
contrast provides a promising technology to meet the increasing performance goals in this 
area. 

Previous experiments done at the Naval Postgraduate School have shown that the 
infrared electric vector of sea radiance near grazing angle appears perpendicularly 
(horizontally) polarized within regions of direct sun glint and parallel (vertically) polarized 
otherwise [Ref. 1,2,3]. Sky background and manufactured targets appear much less 
polarized than sea background. Therefore, polarization filters aligned horizontally while 
viewing target scenes outside the sun glint can improve the target-to-background contrast. 
In addition to these works, concerns about the time dependency of contrast results led to 
the design of a new device to split the scene in two different polarization fields in the same 
image [Ref. 4]. 

In this project data from an experimental data set obtained in the Electro Optic 
Propagation Assessment in the Coastal Environment (EOPACE) multinational 
measurement series in San Diego Bay in 1996 is analyzed. EOPACE is a program 
sponsored by ONR and organized by SPAWAR SYSTEMS CENTER, San Diego, Ocean 


and Atmosphere Sciences Division, Code 543. The overall purpose of the EOPACE 


l 


measurements is to quantify infrared (IR) propagation characteristics for near ocean 
surface transmission and analyze electro-optic (EO) systems performance in the coastal 
environment condition'. The basic data analyzed in this thesis were collected by the NPS- 
Meteorological Group on R/V POINT SUR, NPS-Naval Academic Center for Infrared 
Technology (NACIT) on R/V POINT SUR and NPS-NACIT at Building 15 of the Naval 
Command Control and Ocean Surveillance Center (NRAD) in April 9 and April 10, 1996. 
The data includes 5,508 selected images of the R/V POINT SUR in sets of different 
ranges and aspect angle, meteorological conditions, Global Positioning System (GPS) data 
and ship skin temperature. The collected image sets are sequences of horizontally 
polarized, vertically polarized and unpolarized images taken using an AGA780 camera 
with internal polarizers. With a greater number of images in comparison to the previous 
work done, the main objective of this project is to analyze the effect of the target aspect 
angle on polarized target-to-background contrast. In addition, previously reported results 
related to target-to-background contrast improvement using a horizontally oriented 
polarizer will be verified. Assumptions used in those projects about time dependency of 
contrast results and about generation of unpolarized images by averaging horizontally and 
vertically polarized images will also be verified. 

Chapter II provides the fundamentals of infrared theory. Chapter III describes the 
relevant conditions to obtain and analyze the data. Results are presented and discussed in 


Chapter IV, and finally, conclusions and recommendations are provided in Chapter V. 


! Detailed information about EOPACE measurement series can be found at the 


internet address www.sunspot.nosc.mil 


II. FUNDAMENTALS 


In order to build a gradual approach to the objectives of this work, presenting the 
basic infrared theory is very important. Naturally, there is no intention to develop the 
topics in full detail. However, the most important points to understand the phenomena 
associated with infrared radiation will be shown. 

А. ELECTROMAGNETIC SPECTRUM 


Energy can propagate by spatial and temporal variations of electric and magnetic 


fields. The relationship between these two fields is represented by Maxwell’s equations: 


ҮхЕ----- 
ді 


ар 
VxH-J ir (2.1) 
V:D=p, 


У-В-0 


Where: 


E = Electric field intensity vector (V /m) 

D = Electric flux density vector ( C / п?) 

H = Magnetic field intensity vector ( A/m ) 

B = Magnetic flux density ( T ) 

t = Time (s) 

p," Volume charge density ( C / m) 

Therefore, “electromagnetic waves" such as optical, infrared, radar and radio 


waves have the same nature, being composed of two component fields: electric and 
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magnetic. Basically, different names are used to distinguish different frequency intervals in 
the electromagnetic “spectrum.” 

The infrared region of the electromagnetic spectrum lies between the visible 
(shorter wavelength) and the radio (longer wavelength) regions. Figure 2.1 shows this 
region of the spectrum. Its bandwidth is approximately from 0.7 um to 1000 um. 


However, for conventional IR imaging the practical bandwidth is approximately from 


1.5ит 10 15ип. 
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Figure 2.1 - Electromagnetic Spectrum [Ref. 10] 
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Based on factors such as the form of radiation, methods of detection, and the 
atmospheric transmission, the infrared part of the spectrum is subdivided into several IR 
bands. The most common bands are the near-infrared (NIR) 0.7um - 0.9um , the short- 
wave infrared (SWIR) 1.0um - 3.0um , the mid-wave infrared (MWIR) 3.0um - 5.0um, 
and the long-wave infrared (LWIR) 8.0um - 12um. 

Detectable radiation from 0.7um - 3.0um is mostly that reflected from the target 
scene. Night Vision Devices use reflected radiation and can operate in the IR band limited 
to 1.5um by the detection technology. Infrared radiation with a wavelength beyond 
3.0um is called thermal radiation because the emission, dependent on the target 
temperature, gives most of the detectable signature. Infrared systems such as FLIR 
(Forward Looking Infrared), IRLS (Infrared Line Scan), and IRST (Infrared Search and 
Track) use thermal radiation commonly in the MWIR or LWIR band. 

B. INFRARED TERMINOLOGY AND UNITS 

Since electromagnetic energy comes from a source through the atmosphere until it 
reaches a sensor, the fundamental radiometric quantities are associated with the 
characteristic spreading of the infrared radiation. Therefore, they are related to energy, 
area, solid angle and wavelength units. Table 2.1 shows fundamental radiometric 


quantities and their units. 


TABLE 2.1 


RADIOMETRIC QUANTITIES 





Radiant Flux Radiant Flux leaving an element | Watts/cm' 
Density of a surface per unit area. 
(Exitance) 


E Irradiance Radiant power per unit area Watts/cm' 
incident on a surface. 
Radiant Radiant power leaving a point Watts/Sr 
Intensity source per unit solid angle. 


Radiance Radiant power leaving or Watts/cm’S 
arriving at a surface at a point in 
a given direction per solid angle 
per area normal to that 
direction. 


M characterizes the power leaving the target surface per unit area. Some of this 







radiation is emitted by the target and some is reflected due to /rradiance from another 
source. The distant transport of the radiant energy from the target 1s described by the 
Radiant Intensity, which includes the geometrical spreading with range. When the 
radiation reaches the detector, the /rradiance at the receiver of the target radiation 1s 
converted into a signal for recording. 
C. LAWS OF THERMAL RADIATION 

1. Васк Воду 

When radiation is incident upon a target, fractions of the total radıant energy are 


absorbed, reflected, and transmitted, so that 


1 =0+р+т (2.2) 


where each coefficient represents the respective absorbed, reflected, and transmitted 
fractions. Kirchoff's law of electromagnetic radiation states that a good absorber is also a 


good emitter of radiation by the equation: 


(-Е (2.3) 


A black body is known simply as an ideal emitter and absorber. Therefore, it has 
absorptivity (a) = emissivity (€) = 1, so the reflectance (p) and transmissivity (7) = 0. For 
a black body this is true at all wavelengths. The spectral radiant emittance curve for a 


black body is shown by Figure 2.2, and is described by Planck's Law. 


Wai 


Wary 
em? um 





Figure 2.2 - Black Body Radiant Emittance [Ref. 7] 
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2; Planck's Law 
Max Planck theorized that atomic oscillators do not emit radiation continuously 
but in discrete amounts of energy (quanta). The frequency of the emitted radiation is 


related to the quantized energy change by the equation: 


E=hv (2.4) 


Where 
E = Energy change (J) 
h = Planck’s constant (h = 6.625 х 103 75) 
v= frequency (Hz) 
From these principles, the equation to represent a black body spectral exitance as 


function of its temperature was derived. This equation is known as Planck's Law 


_ 2пс?Һ/А? 
Erna meer (2.5) 
(е 2 -1) 


where: 
M, = Spectral exitance (W / cm? / um) 
4 = Wavelength ( um) 
k = Boltzmann’s constant (k =1.38 x 107 J/K) 
c = Velocity of light in the vacuum (c = 3 x 10° m/s) 
ћ = Planck’s constant (h = 6.625 x 10% Js) 


T = Temperature of the source (K) 


Figure 2.2 shows the graphical representation of Planck’s law. 


For non black body emitters, the radiant exitance is reduced by the emissivity (€). 


The radiant exitance is then: 


М = МА Бтаск-Боду (2.6) 


Where: 


‚ = spectral emissivity (e < 1 for all non black bodies) at wavelength 4. 


Some examples of emissivities of practical materials are shown in Table 2.2. 


TABLE 2.2 


EMISSIVITIES OF COMMON MATERIALS 













If the spectral emissivity (€,) is a constant, then the body is called a gray body. 


Many materials behave as gray bodies over a limited range of wavelength. 


3. Wien Displacement Law 
Wien's displacement law can be derived by differentiating the radiant exitance M, 
(Eq. 2.5) with respect to wavelength and setting the result equal to zero. This procedure 


will give the wavelength of maximum exitance as a function of the temperature: 


p T=2898. 78 (2.7) 


Where: 
Ama Wavelength for maximum exitance (um) 
T - Temperature of the source (К) 
From Wien's Displacement Law it 1s observed that as the temperature increases, 
the wavelength corresponding to the maximum exitance decreases. 
4. Stefan-Boltzmann Law 
If the radiant exitance is integrated over all wavelengths for a given temperature, 


the total power emitted per unit area can be found. The equation is found to be, for a 


black body: 
DUE e (2.8) 
Where: 
o = Stefan-Boltzmann constant ( ¢ =5.67 x 10? W / cm? K^) 


7 = Temperature of the source (K) 


For a gray body the result must be multiplied by the emissivity. 
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D. ATMOSPHERIC PROPAGATION 

The atmosphere is composed by many different elements. Therefore, it emits 
radiation and interferes with the radiation passing through it. The spectral radiant flux 
leaving a source at distance R from the sensor will be attenuated by the atmosphere. The 
attenuation is a function of the distance R and the wavelength of the electromagnetic 
wave. The atmospheric transmission coefficient t (0 < т < 1) represents the general effect 
of the atmospheric attenuation. t decreases as the distance increases. On the other hand, 
the variation of t with the wavelength is a more complex function and depends on the 
atmospheric composition. Several standard models are used to predict the influence of the 
atmosphere. LOWTRAN and MODTRAN are examples of those models. Figure 2.3 
shows a characteristic graphic result for atmospheric transmission coefficient generated by 
MODTRAN for a Mid-Latitude Summer atmosphere, 10 km horizontal path at 100 m 


altitude over the open sea. 


NTNMODSPHERIC TRANSMISSION 
tooo 


15000 3000 3000 2000 


Eccli 
Trane TT 





> 


wvavsliongth Micronz 


Figure 2.3 - Atmospheric Transmission 
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Figure 2.3 shows two intervals clearly better in terms of atmospheric transmission: 
3.0um - 5.0um (MWIR) and 8.0um - 12um (LWIR). 

The attenuation dependency on wavelength is an intrinsic characteristic of the 
atmosphere components associated with two basic mechanisms: scattering and absorption. 

The scattering loss mechanism is the redistribution of incident radiation over all 
directions. The predominant species in this process are suspended aerosol particles for 
which the particle size is close to the wavelength. Aerosols include dust, fog, salt, ash, and 
many other airborne particles. Particularly in the marine boundary layer, aerosols include 
condensed water droplets. The principal effects of scattering are observed at the shorter 
wavelengths in the IR spectrum (below 1.2um). Above 1.2um molecular absorption 
dominates. 

Absorption mechanisms are due to the structure of the molecules and the modes in 
which they can be excited and thus absorb energy. The predominant absorbing species in 
the atmosphere are molecules such as H,0 and CO,. Basically, they define the observed 
"windows" used in IR systems (MWIR and L WIR) as seen in Figure 2.3. 

E. POLARIZATION PHYSICS 

Electromagnetic waves are composed of two fields: electric and magnetic. For a 
plane wave propagating in the z-direction, the electric field components in the plane 
perpendicular to the direction of propagation may be represented by: 


Е (2,0) =E,,cos(wt +6) 


E (2,1) =E,,cos(wt+0,) (2.9) 


1 


Where: 

Е, = x-component amplitude 

E,,= y-component amplitude 

д, = phase constant of x-component 

0,= phase constant of y-component 

The plane containing the normal to a boundary surface and the direction vector of 
a propagating wave is called the “plane of incidence.” The direction of the electric field 
vector in the plane perpendicular to the propagation defines the type of polarization. If the 
electric field vector is parallel to the plane of incidence, the polarization is called 
“parallel.” Similarly, if the electric field vector is perpendicular to the plane of incidence, 
the polarization is called “perpendicular.” In this project, considering the sea surface as a 
flat horizontal reference, parallel polarization related to it will be called “vertical 
polarization.” Similarly, perpendicular polarization will be called “horizontal polarization.” 

Polarization occurs in the reflected and emitted IR radiation of targets and 
backgrounds. Reflection and refraction of light are described by the Fresnel equations and 


Snell’s law. The incident, reflected, and transmitted waves are shown in Figure 2.4. 
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Figure 2.4 - Wave Incidence Geometry 


ІНЕ 


For plane waves incident at a boundary between two media of different complex 
refractive indices n, and n, , the angle of reflection equals the angle of incidence. The 


angle of transmission @, as a function of the incidence angle @, is given by Snell's law: 


n,sin(0) -7,sin(0)) (2.10) 


For nonmagnetic media, the Fresnel equations can be written: 
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The degree of polarization of the target (or background) can be calculated by: 


<N >-<N,> 
Mole (2.12) 


<N > +<N > 


Where: 


<М > = mean radiance with vertical E-field 
<N,> = mean radiance with horizontal E-field 


For the purposes of this project, three different elements are considered as 
components of the infrared images: ship (target), sea and sky (background). Each of them 


has its own polarization characteristics in emitted and reflected radiation. 
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1. Sea Emission and Reflection 

Sea surface polarization characteristics have complex mathematical models. 
Extensive studies have been developed to model the sea wave dynamics to improve 
predictions about the polarization behavior of sea radiance. Many variables influence this 
process and a probabilistic approach is needed to take into account all their effects. 
References 6, 11 and 12 present detailed description of some models frequently used for 
analysis of the polarization characteristics of the sea surface. Previous experiments done 
by the Naval Postgraduate School have shown that the electric vector of sea radiance 
appears horizontally polarized within the sun glint regions and vertically polarized 
otherwise [Ref. 1,2,3]. Figure 2.5 shows the predicted degree of polarization for MWIR 
and LWIR in and near the sun glint region as a function of solar zenith angle for an 
observer looking in the vertical plane of the sun [Ref. 3]. The negative values of 


polarization outside of the sun glint regions indicate parallel polarization. 


Degree of Polarization 
5 
а 
о 


e 
8 


Degree of Polarization 


& 
8 





-0.040 


170 175 180 185 190 
Observer Azimuth Angle 
Wind Speed 3 M/S Observer Zenith Angle 89 Degree 





MWIR LWIR 


Figure 2.5 - Theoretical polarization in and near sun glint [Ref. 3] 
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21 Sky Emission and Reflection 

The IR radiation from the sky is essentially unpolarized [Ref. 1,2,4]. Therefore, 
use of static polarized filtering components will not contribute to the contrast 
improvement measurements when the sky is viewed directly. The sky emission in the 
image will be subject to attenuation by the polarizing filter to the degree of the unpolarized 
transmissivity of the filter. Figure 2.6 shows the sky radiance as a function of zenith angle. 

Sky radiance increases with zenith angle toward the horizon to that of a black- 
body with atmospheric temperature. This is due to the maximum thickness of the 
atmosphere encountered at that angle. Thus, even the small emissivities of the absorbers 


combine. 
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de Target Emission and Reflection 

Although primary contrast improvement in the sea scenario will be gained by 
filtering of sea emission polarization, the polarization due to emission from the target itself 
must be considered. Painted surfaces such as found on all types of vehicles and ships 
display emission polarization [Ref 14]. Figure 2.7 shows 21 measurements taken on 
various paint samples ranging from very smooth texture through sand-paint mixtures. The 
line for an ideal surface in Figure 2.7 denotes polarization that would be seen from a 
perfect specular reflector. The rougher samples show less polarization than the smoother. 
The amount of polarization depends on the index of refraction of the paint and any 


degradation due to surface roughness. 
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Figure 2.7 [Ref. 14] 


ШЕТІ 


The characteristics of target reflectance follow the same behavior as the sea wave 
facet reflectance. The difference in the magnitude of the polarized components is a 
function of the index of refraction of the surface material ofthe target. The net 
polarization will be dependent on the balance between thermal emission and reflection of 
incident radiance. 
F. TARGET / BACKGROUND CONTRAST 

Depending on application, infrared systems must be able to detect, recognize, and 
identify a target against a background. The difference between the recorded irradiance 
from the target and from the background is essential to accomplish those tasks. 

The numerical measurement of this difference is called target / background 
"contrast." Several methods such as signal processing and polarization filtering can be 
used to improve the target background contrast. 


The absolute contrast (irradiance difference) for any polarization can be written as: 


АМ-<М, МЕ ка” (2.13) 


target 


Thus, the improvement factor for contrast for horizontal filters over vertical filters 
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Relative contrast to the background is most commonly defined as: 
Ne, 
CE e 4 (2.15) 


Noked 
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Contrast improvement then is the ratio between polarized contrast and unpolarized 


contrast, 





C pol 
(om = (2.16) 


ипр 


Previous work [Ref. 3] has shown that target / background contrast outside the 
sun glint corridor can be improved by using horizontal polarization filters. The filters 
reduce the predominantly vertically polarized sea background radiance and also, but to a 
lesser extent, the target radiance. Therefore, an overall improvement in contrast could be 
observed. 

The variation of the target aspect angle was not considered in previous 
measurements of ship-sea contrast. As presented in section E.3 in this chapter, a variation 
in the target apparent radiance is expected to be a function of its aspect angle. In this 
work this factor will be studied using the methods presented in chapter III. In addition, 
conclusions and assumptions of previous work in respect to the overall contrast 
improvement, contrast variation and generation of unpolarized images will be verified 
using a much greater number of analyzed images, and including direct measurement of 


unpolarized images for comparison. 
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IIl. EXPERIMENT DESCRIPTION 


In this chapter, all experimental conditions for data acquisition and analysis will be 
described. Equipment, experiment configuration, data structure, software, and analysis 
programs comprise the experimental conditions for the experiment. The EOPACE data 
acquisition of this study was done in April 1996 in San Diego Bay. The analysis covers a 
total of 5,508 infrared images taken of the R/V POINT SUR in April 9 and April 10, using 
an AGA-780 Thermovision camera mounted at Point Loma. Polarizer filters were used to 
get sets of images in three different polarization cases: horizontal, vertical and 
unpolarized. All images were taken in the LWIR band. Meteorological data, GPS data and 
ship “skin” temperature data were recorded by instrumentation installed in the ship. 
Appendix A shows a map of the data acquisition area (San Diego Bay), the GPS ship 
trajectory and the planned ship trajectory. 

А. ЕОШРМЕКТ 

|: АСА-780 Тћегтоу топ 

The AGA-780 Thermovision is a two-channel infrared scanning single detector 
system. It can operate in both the MWIR and LWIR bands using liquid nitrogen cooled 
InSb and HgCdTe detectors, respectively. The nominal working temperature is 77K. 
Scanning is accomplished by two rotating prisms in each channel which scan the image 
across the detector yielding 1.1 mrad geometric resolution. The input lenses for the long 
wave infrared (LWIR) used for the experiments consist of anti-reflection coated 


germanium. Its optical system is f/1.87, with a 7x7 degree field of view. 
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2: Polarizer Filters 

Internal filtering of the radiation after it passes through the imaging system input 
lenses can be accomplished with the AGA-780 Thermovision Scanning Imaging system by 
using an existing filter wheel. The original purpose of the filter wheel was to hold various 
accessory spectral filters in the optical path before the radiation reached the detector. In 
this experiment polarization filters were installed in the filter wheel. Figure 3.1 shows the 
internal layout of the AGA-780. One of seven positions of the wheel is selected by 
rotating the knob on the front face of the unit. The polarization filters used, manufactured 
by Graseby-Specac; Suffolk, England, each consisted of an aluminum grid superimposed 
on a 9.5mm diameter x 3mm thick round KRS-5 substrate. The performance 
specifications provided with the filters are nearly identical. Appendix B shows the 
performance curves for the two filters installed in filter positions 4 and 5, oriented with 


passage axis vertical and horizontal, respectively. 
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Figure 3.1 Internal Layout of the AGA Thermovision (model with only one 
IR Band). [Ref 8] 
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31 R/V POINT SUR 

The R/V POINT SUR, owned by the National Science Foundation (NSF), is 
operated for the Central California Oceanographic Cooperative (CENCAL) by Moss 
Landing Marine Laboratories, which is on Monterey Bay. It is equipped with the 
necessary navigational, laboratory and mechanical facilities that support biological, 
geological, chemical and physical oceanographic research. Setting and recovering current 
meters, chemical sampling, diving, biological trawling, gear testing and seismic reflection 
profiling are among a number of operations within the ship's capability. In this project it 
was used as a cooperative target for infrared imaging as well as for atmospheric sampling 
functions in the EOPACE experiment. Appendix C presents the ship's general 
specifications and dimensions (R/V POINT SUR Cruise Planning Manual). 

4. Software 

a. CEDIP PTRWIN 
The analog data provided by the scanning detector is digitized for display, 

recording and measurement by PTRWIN version 3.16 from CEDIP. This data acquisition 
and analysis software uses a 12 bit A/D converter for converting the detector signal 
intensity into a real-time image on a computer screen and for digital image recording. In 
addition, the PTRWIN software was used as a visualization tool to select images to be 
analyzed. Although PTRWIN also offers good analysis tools, it works interactively 
analyzing only one image a time. Thus, this software package is not adequate to be used as 


analysis software in this project due the large number of images to be analyzed. 
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b. Interactive Data Language (IDL) 

IDL is a complete computing environment for the interactive analysis and 
visualization of data, integrating an array-onented language with mathematical analysis 
and graphical display techniques. All image analyses in this project were performed using 
programs written in IDL 5.0.3 from Research Systems Inc., Boulder, CO. 

с. ОМТАК РСМОРТКАМ 

PCMODTRAN 3 from Ontar Inc. was used in this project to calculate the 
average atmospheric transmission from the available meteorological data. It can operate in 
MODTRAN 3 or LOWTRAN 7 modes. All calculations were performed in MODTRAN 3 
mode using meteorological data (METOC) collected during the experiment. Appendix E.2 
presents the METOC information for each data analysis file. The following input 
parameters were used in the atmospheric transmission calculations: 

1. Atmosphere model: Midlatitude summer 

2. Aerosol model: Navy maritime 

3. Air mass parameter: 3 

4. Geometry: horizontal path at 10 m above sea level. 

5. Distance between sensor and target: from GPS information 
6. Wind speed: from METOC information 

7. Pressure: from METOC information 

8. Temperature: from METOC information 


9. Relative humidity: from METOC information 
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В. EXPERIMENT CONFIGURATION 

The experiment was set in San Diego Bay at Point Loma. Appendix A.1 shows a 
map of the area. The basic components in the experiment configuration were the sensor 
(camera AGA 780) and the target (R/V POINT SUR). The sensor was located at building 
15 in NRaD, Point Loma. Its geographic coordinates were N 32° 39' 36", W 117° 14' 22". 
The ship positioning was planned to correspond to a pattern of radial bearings with 
stationary points spaced at 0.5 nm intervals. At each stationary point, the ship should 
maneuver to be at different heading for each image sequence (3 polarization cases) taken. 
Both ship GPS location and sensor position are known, so the bearing is computed. From 
ship heading and bearing information, the target aspect angle can be calculated. Appendix 
A presents the planned ship positioning for the experiment. 
C. ORIGINAL DATA 

1: Images 

The sets of images were recorded sequentially for three different polarization 
cases: horizontal, vertical and unpolarized. Each polarization case has an average of 26 
images (frames). Horizontal and vertical cases were recorded in the same file whereas 
each unpolarized case was recorded in a separate file. The frames are composed of 4 
interlaced fields with dimensions 125x64 pixels each. Therefore, each frame is an array of 
125x256 pixels. A total of 5,508 images were analyzed. Image data taking was 
implemented with the AGA 780 mounted on a remotely controlled pan-tilt head outside 
the south looking window of Building 15 in Point Loma. The data were digitalized and 


recorded with the PTRWIN acquisition boards installed in the docking station of an NEC 
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laptop computer. For storage, the files were transferred to a second portable computer 
and stored on optical disk. The original image files from PTRWIN are identified by the 


extension *.PTW. Table 3.1 presents the *.PTW binary file structure. 


TABLE 3.1 
ORIGINAL IMAGE FILES (*.PTW) STRUCTURE 
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Header 
(1,536 bytes) 












Comments 





Ist Field 









2nd Frame 
(64,000 bytes) 





2nd Field 16,000 








3rd Field 32,000 





4th Field 48,000 





The Frame Structure above is repeated until the last frame in the image sequence. 


Appendix D.1 presents a list of the selected *. PTW files. 
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2. Meteorological Data 

Meteorological data were collected by the NPS Boundary Layer Meteorology 
Group on 2 meteorological towers (METOC 1 and METOC 2) installed on the ship. The 
data were recorded at 30 second intervals. Ship GPS data are included in the METOC 
files. Appendix E.1 describes and presents the format used to record the data. METOC 1 
tower was located above the ship bridge (38 ft above sea surface) and METOC 2 tower 
was located at the ship bow (24 ft above sea surface). Meteorological conditions for each 
analysis file are presented in Appendix E.2 whereas timing and positioning data are 
presented in Appendix E.3. 

3. Ship Skin Temperature Data 

The temperature of the ship "skin" was measured by 16 thermistors installed by 
NPS-NACIT in different points of the ship. The data were recorded at a 20 second rate. 
Appendix F presents the position for each thermistor and the format used to record the 
data. 
D. PRE-PROCESSED DATA 

All files recorded under the PTRWIN program were stored in *.PTW format. 
The original format used to store the images (*.PTW) has some inconveniences for use as 
a data base for the analysis programs: 

- Polarization cases for one test condition are divided in two different files. 

Consequences : programs would have to open two different files for analysis, increasing 
the time and memory resources needed. 


- Headers have much information that is not used for the analysis. 


20) 


Consequences : increase in readıng time and memory resources needed. 

- Field interlace sequence is not constant for all files. 
Consequences : before each analysis the interlace sequence should be verified and 
corrected by the user interactively for each analysis phase, increasing the analysis time. 

- File structure is specific for the camera settings. 
Consequences: the analysis programs could be used only for this specific format. 
Converting another format (a single JPEG or BMP image for example) to the original 
image format (*.PTW) is not easy. 

In order to solve these problems, all original image files (*.PTW) were converted 

into a new format (*.PTE) with the following properties: 

- All polarization cases included in one file (each file represents one test 
condition). 

- Images are integer arrays (in this project their dimensions are 125x256). 

- Only the essential information for the analysis programs is included. This 
information is given as an integer matrix. 

- All analysis results are included in the same data file as a floating point 
matrix. 

- Converting any image format into a ^*.PTE" file format is easy. 
Therefore, the analysis programs do not need modification to be used to analyze images 
recorded in other formats. 


The file structure for “* PTE” files is described in the next item. 
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1. Images 

140 selected original image files (*.PTW) were converted into 70 "basic image 
files" (*. PTE) using the program BASEGEN.PRO (appendix H) written in IDL. 
Appendix D.2 presents a list of the generated basic files. The denomination “basic file” is 
used because they are the “data base” for the analysis programs. The basic files (*.PTE) 
were named sequentially using the root name “BASE” followed by a number. For 
example, the basic file corresponding to the 10" analyzed test condition (3 polarization 
cases) is BASEIO.PTE. This technique was adopted to simplify the reading of all files 
sequentially. All analysis programs will recognize basic file structures for any dimensions 
of image arrays and any number of polarization cases. 

The basic files structure for 3 polarization cases is presented in Table 3.2. 


TABLE 3.2 
BASIC FILES (*.PTE) STRUCTURE 














! In IDL, the first matrix index is column number and the second is line number. 
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Table 3.3 presents the structure for the information vector variable referred to ın 


Table 3.2. 


TABLE 3.3 
INFORMATION VECTOR 
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Time matrices have 4 columns to describe the time for each frame recording: hour, 
minute, second and hundredths of second.. 

2. Meteorological Data 

The original meteorological data were converted into a float matrix with 
dımensions 9 x 2882 x 2 where the 9 columns represent the average values for 


measurements performed on the 2 different meteorological towers on the ship. Each line 
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corresponds to a data sample and each layer corresponds to one day (April 9 and April 
10). Table 3.4 presents the parameter's columns in this matrix ((METDATA.PTE"). 


TABLE 3.4 
COLUMNS OF THE METEOROLOGICAL DATA MATRIX (METDATA.PTE) 
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Air pressure (mb) 1018.35 
Sea surface temperature (C) 17:30 


Another matrix (“‘AVEMET.DAT”) was also created to summarize meteorological 





data, ship positioning and average atmospheric transmission for all 70 basic files. These 
matrix data are presented in Appendix E.2. Table 3.5 presents a sample of the 


AVEMET.DAT data (BASE01.PTE). 


S 


TABLE 3.5 
COLUMNS OF THE AVEMET.DAT MATRIX (SAMPLE BASE01.PTE) 


ООО 





File number 
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12 Theta (angle between the vertical vector fixed at the sensor and 1.62 
the sensor to the ship vector in radians) 
13 241 


Phi (angle between the True North pointing vector fixed at the 
















sensor and the horizontal projection of the sensor to the ship vector 
in radians) 


Atmospheric Transmittance 





3. Ship Skin Temperature Data 

The original skin temperature data recorded by NPS-NACIT on the R/V POINT 
SUR were converted into a floating point matrix (“SKINDATA.PTE”) with dimensions 
17 x 4321 x 2. Its first column represents the elapsed time in seconds from 00:00 h local 
of the corresponding day. The 16 other columns represent the temperature measured by 
each thermistor. Each line corresponds to a data sample and each layer corresponds to one 
day (April 9 and April 10). 
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4. Thermal Image Calibration 

During the image acquisition process, the detector signal related to a given point 
of the scene is converted to a digital number at the correspondent point in the image. The 
digital number recorded by the digitalization software is called “thermal level.” Although 
the intensity of radiation received at the detector is presumed to be proportional to the 
system's thermal level, this signal is a nonlinear function of temperature. Thus, in order to 
have a correct correspondence between thermal level and temperature, the AGA 780 
requires calibration with a laboratory blackbody source. The system output of thermal 
level is recorded against the source blackbody temperature through a range of 
temperatures and the calibration curve is generated through a least squares fit method and 


is of the form [Ref. 8]: 


TE OS (3.1) 


Where: 

- A, B, and C = Constants determined by the curve fitting 

- OS = Offset correction to adjust for detector response drift [Ref. 1] 

- T = Temperature (K). 

Appendix G presents the laboratory calibration curves for the 3 analyzed 
polarization cases. Thermal level values are used to calculate the image contrast. The 
camera was set to atmospheric transmittance and target emissivity equal to 1.0. Sınce the 


images are taken at different distances and the target is not a blackbody, their thermal level 
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values must be calibrated before the analysıs process. Since the skin temperature was 
available, the image calibration procedure described in the AGA Operating Manual using a 
known temperature reference source was used. In this process, the calibrated thermal level 


is calculated by the following equation: 


АГ, 
Г. = i (3.2) 





Where: 
- AL, - Difference between a point thermal level and the recorded thermal level at 


the reference point 


-T, = Average atmospheric transmission 
-6 = Target Emissivity 
- L — Calibrated thermal level at the reference point in the target (eq. 3.1). 


The calibration process can be described as follows: 

- The recorded image 15 composed by a matrix of uncalibrated thermal level 
values (one for each image pixel). 

- The stored thermal level values are assumed to be proportional to the 
apparent radiance at each given pixel. 

- The skin temperature at a reference point (stack) on the target is known. 

- From the laboratory calibration curve for that polarization case, a thermal 
level value corresponding to the reference point temperature is calculated and the 


calibrated thermal level for any point is calculated using Equation 3.2. The thermal level 
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value at the reference point, calibrated by the calibration curve, is called “reference 
thermal level (L)." It considers the target as a blackbody and it is related to the real 
temperature (it does not consider the distance between the sensor and the target). The 
idea behind the calibration process 15 to relate each thermal level value in the image to the 
reference thermal level conditions (blackbody at short distance). Although the 
mathematical derivation of the Equation 3.2 can be followed in the reference 8, its physical 
meaning can be explained using simple ideas: 

- The calibrated thermal level at the reference point (reference thermal 
level) considers the target as a blackbody at short distance. 

- To be calibrated, the recorded thermal level at any point is going to be 
brought at the same conditions of the reference thermal level (blackbody at short 
distance). 

- The target is assumed to be a gray body. Therefore, to correct the 
recorded thermal level of the target points for the blackbody condition, it is necessary to 
divide its values by the target emissivity. 

- Since the recorded thermal level is assumed to be proportional to the 
radiance, to correct it for the target distance to the sensor, it is necessary to divide its 
values by the atmospheric transmission coefficient. 

- Since the reference point is in the target, the path radiance has the same 
value for the reference and the other points. Therefore, by subtracting the reference 
thermal level from a point calibrated thermal value, the path radiance is canceled and thus, 


it is not present in Equation 3.2. 
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Е. ANALYSIS PROGRAMS 
The purpose of the analysis carried out was to analyze statistically the sequence of 
frames with each polarization for a given case, so as to compute; 
- Averaged degree of polarization for ship, sea background and sky 
background. 
- Temporal variability of each and temporal trend if any. 
- Target to background contrast for each polarization. 
- Contrast improvement for horizontal and vertical polarization cases. 
- Degree of polarization of the ship and of its major planes as functions of 
aspect angle. 
To carry out this analysis it was necessary to develop a set of computer programs. 
All programs were written in IDL. The complete set comprises 37 different programs that 
can be divided in 3 categories: 
- Auxiliary programs: specialized functions used to perform common tasks 
for other programs. 
- Data management programs: programs used to perform data format 
conversion and results summarization 
- Executive programs: programs used to calculate results from the available 
data. 
Appendix H.1 presents a summary table of the program's name, type and 


objective. 
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E, Methodology 

All programs are completely commented in their text sources to help the 
user to understand and modify them for future analysis’. The set of programs runs in a 
sequence given by PONTES.PRO program for each image sequence (polarization case): 

- HOTSPOT.PRO: Finds and Tracks the ship hotspot (stack). 

- AREAGEN.PRO: Defines an “analysis box” around the ship. 

- HORIZON.PRO: Determines the horizon line equation. 

- NORMALIZE.PRO: Normalizes the image to threshold the ship. 

- ELEMENTS.PRO: Isolates each image main component (ship, sea, sky). 

- ASPECT.PRO: Estimates the ship aspect angle. 

- CONTRAST.PRO: Calculates ship (and its major planes) to background 
contrast and degree of polarization for each image main components. 

Although they are commented in their text sources, because of special 
characteristics or methodology, it is interesting to make complementary comments about 
some of them. 

а. Hotspot.pro 

HOTSPOT.PRO is the first executive program in the analysıs sequence. It 
has the basic interactive concept used for all programs which need any kind of user initial 
input. In this concept, the user will be asked to provide initial input data for the first frame 


of the set. The program performs the calculation for all other frames of the set and then 


* All programs used in this work are available with Prof. Alfred W. Cooper, PhD 


at NPS Physics Department 
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shows its results graphically. The user can either accept the computer results or modify 
them (any frame). Finally, the program shows the final result, including the user’s 
modifications. For HOTSPOT.PRO, the initial data requested from the user are the 
definition of a rectangular area using the mouse on the first frame image. The program 
searches for the hottest point inside that area and records its thermal value and position. 
For the next frame it translates the center of the search area to the position of the previous 
hotspot (last frame). This technique allows usage of cluttered images (other “hot” targets 
in the image) for analysis and permits the program to keep the “tracking” on the ship 
hotspot (stack normally) even if it is in motion. Appendix H.2 presents a sequence of 
frames showing a hotspot in a moving target being tracked in a cluttered image. 

b. Areagen. pro 

AREAGEN.PRO generates the coordinates for two corners (diagonal end 
points) of a rectangle around the ship. This rectangle is called the “analysis box” and will 
be used as image object by the other analysis programs. This technique allows the analysis 
to be performed only in the neighborhood of the ship, excluding any possible secondary 
target and so permitting cluttered images to be used. The presence of a secondary target 
inside the analysis region can cause poor results because the thresholded "target" pixels 
might include the secondary target, altering the real target characteristics. Alternatively, 
the secondary target could be computed as “background” pixels, giving a false higher 
average background radiance. AREAGEN.PRO uses the hotspot position as a reference to 
locate the analysis box for each frame. Therefore, if the ship moves, the box will move 
with it. Appendix H.3 presents a sequence of frames showing the box moving along a 


target in a cluttered image. 
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E Horizont. pro 

HORIZONT.PRO calculates the coefficients for the horizon line equation 
by linear regression using 3 points at the horizon line: one point at the middle “x” value of 
the image and two at extreme “x” values (0 and x size of the image). For each point, the 
“y” position is calculated using the natural change in the image thermal level gradient in 
“y” direction at the horizon plane. This change in the gradient can be seen as a “wall” ın 
the surface plot of a typical image in appendices H.5 and H.6. Immediately above the 
horizon, the sky radiance is much higher than the sea radiance. This difference in radiance 
is detected by the program during the process of horizon “y” position searching. The 
initial user input is used to limit the searching interval, assuring that secondary targets at 
the edges of the image will not interfere with the gradient verification process. The idea is 
similar to the hotspot searching, but unidimensional. 

d. Threshold.pro 

THRESHOLD.PRO is used to isolate the ship from the background. The 
basic concept used in this program is the sudden increase in the number of “target pixels” 
(pixels with thermal value above a "testing threshold level”) when the "testing threshold 
level” drops into the background level. In order to apply this concept, the interval between 
the maximum and the minimum thermal values in the image is divided into 100 “testing 
threshold levels.” Starting from the highest testing threshold level, the number of pixels 
above the threshold value is computed and recorded for each one of the 100 defined 
levels. The number of pixels above threshold as a function of the testing threshold level is 


then analyzed. The biggest positive variation in number of pixels is detected and the 
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optimal threshold level is determined. Although the concept is very simple, some 
important points must be taken into consideration to assure good performance of the 
program in detecting the ship with the minimum shape distortion. The most influential of 
them is about the characteristic high radiation level of the sky. As commented in the 
previous item, the increase in the thermal level immediately above the horizon line is seen 
as a “wall” in the surface plot of a typical image in Appendices H.5 and H.6. The thermal 
level of the sky portion in the image 1s much higher than the thermal level in the sea 
portion. Compared with the ship’s thermal level, the sky’s thermal level is normally 
smaller, but as the distance between ship and camera increases, the sky thermal level can 
become higher than the ship (target) thermal level. This is a clear problem to apply directly 
the threshold concept to the original image. By doing this, it would be very likely to 
reduce the thresholded ship area significantly or even to reduce it to zero (detection of the 
sky instead) when the ship distance was enough to make its thermal level less than the sky 
thermal level. To solve this problem, the program “NORMA.PRO” was created. This 
program uses the natural characteristics of the sky (and sea) thermal level variation in the 
“x” direction (horizontal) to normalize the entire image and to eliminate the difference 
between the thermal level of the sky and sea portion without distorting the shape of the 
ship. Observing the surface plot of a typical image in Appendix H.6, the horizon can be 
identified by the base of the “wall” created by the much higher thermal level immediately 
above the horizon when compared with the sea thermal level. In the same picture, since 
the horizon has been identified, it can be seen that the sky and sea thermal level have 


practically no variation in “x” direction (horizontal). This statement is true except for the 
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image’s lines where for the ship (target) is seen to be. Therefore, the normalızation 
process starts considering the horizontal" direction by dividing each pixel's thermal value 
of the image by the average thermal value of its entire line. This first phase will bring the 
sky thermal level and the sea thermal level near to one for all lines “without target in." 
The target pixels will have thermal level values bigger than one and the background lines 
(sea or sky) “with a target in" will have thermal level values less than one. This effect оп 
the background lines that contain the target is seen as a horizontal “depression” in the 
surface plot for normalization phase one in appendix H.7. In a second normalization 
phase, the same process 1s performed in “y” direction (vertical) dividing each pixel’s 
thermal value of the image by the average thermal value of its entire column. These two 
phases comprise the normalization used during the threshold process. Appendix H.8 
shows the surface plot for normalization phase 2. 

е. Elements.pro 

ELEMENTS.PRO is used to separate the main elements of the image: ship, 
sky and sea. The program uses data from HORIZONT.PRO and THRESHOLD.PRO to 
generate “masks” for each element. “Masks” are images with the same dimensions of their 
object analysis box image, but with pixel values equal to 1 for the corresponding element 
and zero anywhere else. Thus, when a mask image is multiplied by its object image, only 
those pixels covered by the mask (its pixels with value 1) will have the original value. All 
others will be zero. Appendix H.9 shows an example of object analysis box, its masks and 


its corresponding elements separated by the process. 
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T Aspect.pro 

The ship heading information recording system had a problem during the 
data acquisition process so this information was not available in the 30 second rate 
МЕТОС data files. The ship aspect angle ıs necessary data to the analysis performed in 
this project. Thus, the program “ASPECT.PRO” was written to estimate the ship aspect 
angle using user provided initial information to start the estimation process which is based 
on target / model shape correlation and aspect ratio comparison. The initial information 
about the ship's aspect angle comes from the notes taken during the experiment by the 
imager operators by radio contact with the ship. These notes contain information about 
polarizer used, ship magnetic bearing and ship magnetic heading. From this information, 
an estimated ship aspect angle can be calculated. ASPECT.PRO searches for the best 
matching in shape and aspect angle turning a scaled model 20 degrees around the initial 
estimation. Calculated results are then shown to the user to be verified and modified if 
necessary. 

The model used by this program 15 a scaled model of the ship composed of flat 
surfaces (planes). A total of 32 planes defined by 76 vertices comprise the model. The 
vertices are defined by three coordinates (x,y,z) and their values correspond to the ship 
dimensions in feet. Although only the four main planes that define the hull sides of the 
model are used for contrast calculations, the other planes used to model the upper decks 
are necessary for the shape matching process used by ASPECT.PRO. The model is 
generated by “MODELE.PRO” and “DRAW_MDL.PRO” and it can be viewed at any 


combination of elevation and azimuth angle. Appendix H.10 presents the model viewed at 
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different angles. The model can also be drawn as a wıreframe structure. The wireframe 
model is used to superimpose the model on the ship image to verify the shape matching 
and model positioning. Appendix H.11 shows the wireframe model superimposed on a 
ship image. 

2. Contrast.pro 

CONTRAST.PRO is used to calculate the contrast for the ship and its 
“main planes.” The ship’s “main planes” are the four biggest planes that comprise the 
- superstructure of the ship model. Appendix H.12 presents a picture showing the ship’s 
main planes. The ship/background contrast ıs calculated by considering the contrast 
between the ship's image portion above the horizon against the sky and the ship's image 
portion below the horizon against the sea. In order to calculate the contrast for each main 
plane separated, the scaled model in “mask” mode is turned to match the ship aspect 
angle. "Mask" mode is an available feature of “MODELE.PRO” which generates the 
model as a “mask” for a given plane. The “mask” is just a part of the image where the 
pixels have value equal to 1. All other pixels have value equal to zero. Therefore, the 
process to isolate a part of an image using a mask image is described as follows: 

- A ship image is represented by an integer matrix where each element 
represents a pixel. 

- Suppose it is necessary to isolate the main “plane 0" of the ship image. 

- The ship model is turned to match the ship image in size and aspect angle. 

- The ship model in mask mode is an image represented by an integer 


matrix where only the pixels corresponding to the main “plane 0" have value equal to 1. 
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All other pixels (image matrix elements) have value equal to zero. 

- By multiplying the ship image matrix by the ship model matrix (mask), 
only those pixels in the ship image corresponding to the main “plane 0" will keep their 
original values. All other pixels will have value zero. Therefore, by this process the main 
“plane 0" is isolated from the rest of the ship image. 

Appendix H.13 presents a figure showing plane mask being applied on a ship 
image. For each of the 4 main planes, the contrast is calculated against the average 


background. Depending on the ship aspect angle some planes are not visible. 
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IV. RESULTS AND ANALYSES 


In this chapter, all experimental results will be described and analyzed. In addition, 
the IDL programs used in this project will also be analyzed. Although these programs 
were used as an analysis tool in this project, they comprise a very useful package for 
future analysis. Thus, improving their performances and interfaces is important. 

А. IDL PROGRAMS 

The programs have shown good performance in terms of stability and quality of 
results. However, it was observed that the execution time and the user interface still need 
improvements to be used as an analysis package. The average time to process one file 
(*.pte) is 40 minutes and there are too many keyboard inputs requested from the user. 
Therefore the following improvements are recommended for future versions to optimize 
an analysis package: 

- Use of “WIDGETS”’ to integrate the programs. 

- Optimization of array treatment to reduce the number of “loops” in the 
programs. 

- Inclusion of “PTRWIN” (from CEDIP) as an external program that can 


be accessed from the main program in the package. 


' IDL graphic tools. These tools permit generation of “windows” type 


applications. 
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В. RESULTS DESCRIPTION 

The results are systematically stored during the analysis process in three matrices 
(one for each polarization case). These matrices are called "result matrices" and have 100 
columns and a number of rows equal to the number of frames for the analyzed polarization 
case. Appendix I.1 presents the format (column description) of the result matrices. 

After analyzing all 70 files, their corresponding result matrices (three each file) 
were summarized by averaging the frame values. The summarized values were then 
- collected in one summary matrix with 130 columns (parameters) and 70 rows (files). The 
format of this matrix (column description) is presented in Appendix I.2. The summary 
matrix 15 too big to be presented completely. However, its contents are presented partially 
in Appendix I in the tables referred to in the next analysis items. 
С. SHIP CONTRAST VARIATION IN TIME 

In this and in previous projects, image sets of different polarization were taken 
sequentially. This means that the image recording starts with the horizontal polarizer 
setting for approximately 8 seconds. Then, the polarizer is switched manually to the 
vertical setting. Finally the unpolarized case is selected on the camera. This change from 
polarized to unpolarized requires resetting the thermal level range controls and recording 
these in the computer. This process involves a significant time lag as seen in Figure 4.1 


(sample from BASEO1.PTE). 
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Figure 4.1 Example of image set time line (BASE01.PTE) 


Since comparisons between the horizontal and the vertical polarization cases have 
to be made, there were concerns about the validity of these comparisons because of the 
possible time variations in the contrast during recording time for those cases. To verify the 
validity of those comparisons, the contrast variation for the three polarization cases was 
analyzed by calculating the contrast for each frame. The results were plotted against the 
elapsed time and a straight line was fitted to the points for each case. The basic parameters 
to be observed in this analysis are the angular coefficient and the correlation coefficient for 
the line fitting. The angular coefficient shows how the contrast increases or decreases in 
time. Thus, an angular coefficient near zero would represent that the contrast does not 
vary significantly in the analyzed interval of time. The correlation coefficient shows how 
well the contrast variation can be represented by a straight line and so how well it can be 
extrapolated in time. For each set of frames, the contrast's mean and standard deviation 
are also calculated. These parameters are useful to analyze the contrast fluctuation during 
the recording interval of time. Statistical considerations about the results are analyzed by 
calculating results for all 70 files and plotting a histogram for each analysis parameter. 
Appendix 1.3 presents the tabulated results for the curve fitting process to the ship 


contrast variation data. Each table shows the results for all 70 files for each polarization 
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case: horizontal, vertical and unpolarized. Appendix J.1 graphically presents a sample of 
the line fitting process for one file (BASE01.PTE). As an example of this process, Figure 


4.2 presents the line fitting results for the vertically polarized case in BASE01.PTE. 


VERTICAL POLARIZATION 
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Figure 4.2 Example of line fitting process (BASE01.PTE) 


Appendix J.2 presents the statistical distribution (histograms) of the fitting 
coefficients for horizontal, vertical and unpolarized cases. Those histograms summarize 
the results calculated for all 70 files. From the histograms it is possible to obtain Table 4.1 


which represents each parameter as a statistical population with 70 samples. 
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TABLE 4.1 
SHIP CONTRAST VARIATION: HISTOGRAMS RESULTS 


Polarization Polarization 
ње [ше |ом fo — 


Table 4.1 : Statistical distribution of curve fitting parameters for all 70 analyzed files.* 
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Standard Deviation 






Analyzing the results, the following observations can be made: 

- The correlation coefficient for the three polarization cases is greater than 
0.8. Thus, the contrast variation in time for the three polarization cases has linear 
behavior. 

- For the three polarization cases, the angular coefficient has mean value 
very close to zero and small standard deviation. This demonstrates that the contrast does 
not have significant variation in time. 

- The analysis of the standard deviation of the contrast distribution shows 
that the unpolarized case has bigger contrast fluctuation in comparison with the vertical 
and horizontal polarization cases. This can explained by the polarizer effect in reducing 


fluctuation of the radiance components for the horizontal and vertical polarization cases. 


^ Appendix E.2 presents all specific conditions for each one of the 70 analyzed files 
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These analyses support the valıdity ofthe results from previous work by 
comparing horizontally and vertically polarized images taken sequentially. 

D. DEGREE OF POLARIZATION AND CONTRAST IMPROVEMENT 

Previous projects have shown that using a horizontal polarizer it 1s possible to 
improve the target to background contrast for sea / ship targets by up to 15% [Ref 1,2]. In 
this project, using a larger number of images and different meteorological environment, it 
15 possible to verify those results with a good statistical confidence. A total of 5,508 
images were analyzed and their results were summarized in 70 files. Appendix I.4 presents 
the target degree of polarization and the target to background contrast improvement for 
horizontal and vertical polarization cases for each one of the 70 analyzed files. 

Appendix J.3 shows the statistical distribution of the contrast improvement factor 
obtained with a horizontal and vertical polarizing filters. Whereas in the MAPTIP 
measurements off the Dutch coast in the late fall the observed contrast improvement 
ranged only above 1 up to 1.5 with the horizontal filter, the EOPACE San Diego 
measurements analyzed here range from 0.5 to 1.8, with a mean of 1.08. This implies that 
use of the horizontal filter to reject the sea surface polarization gives on the average an 
8% contrast improvement, as opposed to 15% in the MAPTIP work. 

The histograms of degree of polarization for the ship target, the sea and the sky are 
also shown in Appendix J.3. For analyzing those results it must be observed that a positive 
degree of polarization shows a predominance of vertically polarized radiation whereas a 
negative degree of polarization means a predominance of horizontally polarized radiation. 


The sea surface shows rather consistent vertical polarization with mean at 0.1218, while 
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the ship shows a very narrow distribution centered effectively at zero. The sky polarization 
is apparently random with widely spaced outliers an both sides. 
Table 4.2 summarize all the results from the histograms shown in Appendix J.3. 


TABLE 4.2 
DEGREE OF POLARIZATION AND SHIP CONTRAST IMPROVEMENT 


Degree of Polarization 


ha — — pm — 
Ship to Background Contrast 
Improvement Standard Deviation 0.242 


Table 4.2: Statistical distribution of degree of polarization and contrast improvement for all 70 analyzed files.’ 

























Sea: 
Degree of Polarization 








Sky: 
Degree of Polarization 












Horizontal Polarization: 
Ship to Background Contrast 
Improvement 












Vertical Polarization: 





From the analyses, the following observations can be made: 
- Ship and sky degree of polarization have mean value very close to zero. 
Thus, neither horizontal nor vertical polarization is predominant in the radiation from the 
sky or from the ship. 
- Positive mean value for sea degree of polarization demonstrates that sea 
radiation is partially vertically polarized over the testing conditions. 


- Mean value for contrast improvement using horizontally oriented 


? Appendix E.2 presents all specific conditions for each one of the 70 analyzed files 
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polarizer shows a gain of 8% (1.081) in contrast in this polarization case. This value is 
smaller than previous results (15%). A higher average sea surface temperature could be a 
factor in this difference. For the analyzed data set, the average sea surface temperature 
was 17.2 C (Appendix E.2) against 11.5 C in previous reports [Ref. 1]. 

These analyses demonstrate that the results from previous work in terms of 
contrast improvement using a horizontal polarizer are qualitatively valid. 
E. GENERATION OF UNPOLARIZED IMAGES 

In previous projects, unpolarized images were generated by averaging the thermal 
level values for horizontal and vertical polarization cases. In the analysis of the Monterey 
Bay and MAPTIP polarization measurements, the Contrast Improvement Factor (relative 
to the unpolarized case) was calculated by constructing an effective unpolarized image as 
the average of the two polarized components. This process is subject to uncertainties 
caused by the presence of the filter with possible transmission extinction loss and thermal 
emission and reflection from the back of the filter. In this project, three polarization cases 
were recorded and analyzed: horizontal, vertical and unpolarized. Therefore it is possible 
to verify the validity of the previous assumption used for composing unpolarized images. 
Appendix J.4 presents plots of the mean thermal level for the unpolarized images against 
the calculated average value using the horizontal and the vertical polarized images. Those 
plots represent the results for the main components of the image: ship, sea and sky 
respectively. Each plot contains 70 points corresponding to the 70 analyzed files. A linear 
fitting process was done to the points in each plot. An angular coefficient close to 1 (45 


degrees) would show that the average assumption was correct. Histograms showing the 


I2 


difference distributions are also presented in the same graphics. In order to illustrate the 
comparison process by line slope, a sample plot for the Sea Thermal Level, the most 
sensitive case, is shown as Figure 4.3. The plot indicates that the average value slightly 


exceeds the unpolarized case, with slope of 0.839 in place of 1. 
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Figure 4.3 Sea Thermal Level 
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Table 4.3 presents the numerical results from the linear fitting process and from the 
histograms in Appendix 1.4. 


TABLE 4.3 
COMPOSITION OF UNPOLARIZED IMAGES 


MA A a a 


a B 
Average - Unpolarized 
Standard Dev. 42.029 38.274 63.461 


Table 4.3: Linear fitting for unpolarized image generation by averaging vertical and horizontal polarization 
thermal levels for all 70 analyzed files.* 





Analyzing the results, the following observations can be made: 
- For all components the mean difference between average and unpolarized 
value is positive. Therefore, the average value is slightly bigger than the unpolarized value. 
- The biggest difference is found to occur for the sea thermal values. This 
fact can be observed by the bigger mean difference value or by the smaller angular 
coefficient when compared with the ship or the sky values. 
- The angular coefficient for the three main components is close to 1.0. 
These analyses demonstrate that the assumption used in previous work in terms of 
unpolarized image generation by averaging horizontal and vertical polarization cases is 


satisfactory. 


* Appendix E.2 presents all specific conditions for each one of the 70 analyzed files 
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Е. SHIP ASPECT ANGLE EFFECTS ON POLARIZATION 

In this project, considering results from all 70 analyzed files, the sea radiation has 
shown a vertical polarization predominance whereas the sky and the ship radiation have 
shown no polarization predominance. The average ship target image has not been shown 
to exhibit any significant degree of polarization in either of the previous measurements 
series. However, literature references have found plane painted surfaces © exhibit 
considerable emission polarization, strongly dependent on angle. It would be expected 
· thus that the ship target would show polarization signature varying with the orientation of 
each facet of the ship surface. Thus the overall polarization of the ship radiance (combined 
emission and reflection) is likely affected by the ship aspect angle. In order to verify this 
relationship the ship was modeled as a set of facets, or planes, as described on item 3.E.1.f 


on page 41 and in Appendix H10 to H12. Figure 4.4 presents a view of the ship model. 





Figure 4.4 Ship Model View 


vo 


The basic idea is to analyze the variation of the degree of polarization of the ship 
and its main planes with their corresponding aspect angle. A major portion of the 
experimentation with the R/V POINT SUR was devoted to evaluation of degree of 
polarization versus aspect angle. The programs ASPECT.PRO, MODELE.PRO and 
CONTRAST.PRO were developed to evaluate separately for each of the major facets of 
the ship the polarization and the contrast for each measurement case. The data for the 70 
files were analyzed with respect to major planes 0,1,2, and 3 as defined in Appendix H.12, 
and also with respect to range. Planes 0 and 3 constitute the port and starboard bow 
quarter of the ship respectively, while planes 1 and 2 represent the port and starboard 
beam. The radiance from these planes was selected using the “mask” procedure explained 
in item 3.E.1.g. Appendix J.5 presents plots of the degree of polarization of the ship and 
its main planes against the aspect angle for different ranges. A variation of degree of 
polarization mostly at the vicinity of the Brewster angle was expected to be observed. 
However, analyzing the graphics in Appendix J.5 it is observed that the points for any 
range are spread around the “x” axis. These results show no strong evidence of any 
functional relationship that could represent a variation of the degree of polarization with 
the target aspect angle in the testing conditions’. The aspect angle has certainly a partial 
influence on the ship degree of polarization. However, over the experiment conditions 
other factors and testing uncertainties can have influence strong enough to cover any 
possible observation of the aspect angle influence on the degree of polanzation. Basically 


three factors in the data processing contributed to the uncertainties in this project: ship 


> Appendix E.2 presents all specific conditions for each one of the 70 analyzed files 
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modeling using planes, ships’ oscillation due to wave movements and lack of ship heading 
information at an adequate sample rate. The aspect angle of each main plane used to 
model the ship's hull certainly can represent the median aspect angle of the corresponding 
part. However, since the hull's surfaces are generally curved, there is an angular error for 
each point on the hull surface. Considering the ship's geometry and the main planes 
positioning, the maximum error due to this factor occurs at the planes O and 3 and it is 
estimated to be 10 degrees. The ships' oscillation is not followed by the ship model. 
Therefore, the angular variations due this factor are not computed in the data process. The 
ship's heading information was available at 10 minutes sample interval. This was not 
sufficient in terms of precision for the image's aspect angle information since each frame 
takes 0.3 seconds. ASPECT.PRO program was then used to estimate the ship heading 
using aspect ratio and shape correlation. In order to test the accuracy of this program, four 
observed headings were compared with the corresponding estimated headings. Table 4.4 


presents the results of the comparison. 


TABLE 4.4 
ASPECT.PRO RESULTS VERIFICATION 





Table 4.4: Comparison between results from program aspect.pro and recorded ship heading information 
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Results in Table 4.4 show an angular error comparable to the plane modeling 
error. The greatest difference (11 degrees) occurred at the longest range (5 km). This is a 
limitation of the program methodology, since it is based on target image contour which 
loses definition with the distance. In order to observe the aspect angle effect, a more 
specific experiment can be designed to reduce the influence of other factors as well the 
uncertainties of this project. For example, a ship carrying a special target could be used for 
this specific experiment. The target would be a metallic plane with known surface 
: temperature. The temperature should be controlled to be constant and higher than the 
ship's stack temperature. A stabilized gyro platform mounted on the target and a GPS 
would give the target's angular data and the ship's positioning to be recorded at least at 1 
Hz sample rate. 

By performing this experiment the uncertainties would certainly be reduced. 
However, the operational meaning of it would be debatable, since it represents conditions 
very unlikely to be found in a real engagement scenario. From this point of view, the 
results of this project, showing no strong evidence that the aspect angle affects the target 
degree of polarization, can be interpreted as a first indication that the target aspect angle 
should not be considered as a primary factor in designing an IR sensor improved by 
polarizers. Another consideration that must still be done is the fact that all results in this 
project are related to the R/V POINT SUR as a target. Therefore, further experimentation 


and modeling work for different targets will be required for final results validation. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The main objectives of this project were to determine a possible influence of the 
target aspect angle on the degree of polarization of the target radiation and also to verify 
previous work results and assumptions related to the target to background contrast 
improvement using horizontal polarizers on the sensor, the time dependency of contrast 
results and the methodology for unpolarized image generation by averaging vertically and 
horizontally polarized images. From the analyses done, the following conclusions were 
drawn: 

1. Although it was theoretically expected, the target degree of polarization 
has shown no strong evidence of any observable dependency of the target aspect angle. 
Therefore, it was concluded that within the experimental conditions other factors and 
testing uncertainties had influence strong enough to cover any possible observation of the 
aspect angle influence on the degree of polarization. In order to observe the aspect angle 
effect, a more specific experiment is proposed, but the operational meaning of it would be 
debatable, since it represents conditions very unlikely to be found in a real engagement 
scenario. From this point of view, the results of this project show no strong evidence that 
the aspect angle affects the target degree of polarization and can be interpreted as a first 
indication that the target aspect angle should not be considered as a primary factor in 
designing an IR sensor improved by polarizers. 

2. The average contrast improvement achieved by using a horizontal 
polarizer on the sensor was found to be 8%. This value is smaller than previous results 
(15%). A higher average sea surface temperature is presented as an important factor for 


explaining this difference. For the analyzed data set, the average sea surface temperature 
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was 17.2C against 11.5 C in previous reports [Ref. 1]. Therefore, the analyses support 
the results from previous work in terms of contrast improvement using a horizontal 
polarizer. 

3. The contrast time variation was not shown to be significant during the 
image recording time. This demonstrated that the results from this and previous work by 
comparing horizontally and vertically polarized images taken sequentially are valid. In 
addition, the analysis of the standard deviation of the contrast distribution has shown that 
the unpolarized case has bigger contrast fluctuation in comparison with the vertical and 
horizontal polarization cases. This fact was interpreted by the polarizer effect in reducing 
fluctuation of the radiation components for the horizontal and vertical polarization cases. 

4. Generation of unpolarized images by averaging horizontal and vertical 
polarization cases has shown to be satisfactory. This was an assumption used in previous 
work. Therefore, the results of previous work are validated. 

Besides the data analysis results, the set of programs used for the analyses have 
been shown to be useful as an analysis package for future work. However, in order to 
improve their performance and user interface, some recommendations were made: 

- Use of “WIDGETS”” to integrate the programs. 

- Optimization of array treatment to reduce the number of “loops” in the 
programs. 

- Inclusion of “PTRWIN” (from CEDIP) as an external program that can 


be accessed from the main program in the package. 


! IDL graphic tools. These tools permit generation of “windows” type 


applications. 
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APPENDIX A 


EXPERIMENT LOCATION 


Description : 


Map ofthe San Diego Bay. 

Plots of the GPS information of the ship trajectory during the experiment in April 9 and 
April 10, 1996. Image set recording points are marked as diamonds in this figure. 

Map showing the planned ship positioning and desired ship heading in each station 


location. 


Conditions : 


The camera (AGA 780) was located at the Building 15 of the NCCOSC-NRaD in Point 
Loma. Its geographic coordinates were N 32 39 35.8, W 117 14 22.2. It was mounted on a table 
projecting from the window, and directed using a remote pan-tilt head. The ship maneuvered in 
reference to the camera position according to five planned magnetic bearings (camera to ship) 
ranging from 135 to 200 degrees. For each bearing there are 5 station locations at ranges from 0.5 
to 3.0 nm. In each station location the ship is supposed to turn to the indicated magnetic headings 


and stay still during the time needed for the image recording process. 


Notes : 


Both ship GPS position and camera location are known, so bearing is computed. Thus, 
ship heading information is used with the bearing information to give the aspect angle of the 
target. 

Range to the shore station and ship heading were reported from the ship at irregular 
intervals and recorded at the sensor. Magnetic bearings to the ship were recorded on land in the 


sensor log. 
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APPENDIX B 


POLARIZERS PERFORMANCE CURVES 


Description : 
Two figures showing the performance curves for the KRS-5 , 9.5mm polarizers used in 


the camera setting positions 4 and 5 respectively. 


Conditions : 

The data were provided by the manufacturer of the polarizers, Graseby-Specac. The 
measurements were taken with a spectrometer using the actual filter shipped. Where two filters 
were required for the crossed-grid measurement (trace 3), another polarizer nearly identical to the 
one shipped is used. The filter performance specifications require 85% transmission efficiency 


over the range 1-12 Am and 98% polarization effectiveness. 


Notes : 
Trace 1 : Transmittance with grid perpendicular to the E vector of the incident radiation 
Trace 2 : Transmittance with grid parallel to the E vector of the incident radiation 
Trace 3 : Transmittance of two filters with grids crossed 
Vertical axis : Wavenumber/wavelength 


Horizontal axis: Transmittance 


65 


АМ TW (d Cis ГГ : Да а 08 Др e de | Ex EE A E. O. AO 
тулоти | BACK OFF С С, — 
ОГ ет 





AS A 
. 


Мə 


ен о > -= , 


ET 
= 


Шер 








ul 








DR 
о ұтты» алаш” 
оа n 


Ма орва ри i] ; 
een а Г. | 15 Aat 
mi E ШҮП {Нит ШЕ li А 
iili: us Malt НЕ 
таға [ETS ПА i [eer Hee 
Eid aee Hil не ЕИ al 
ilr iste an. Е Е 
+ | Ben ИМЕ 
T 4 N ing by! 
... a! A 1 

bo 


- 
E = 
e.c T аво о 


———— фо 
.. 9» 
yo ven 

т. 
be 


е ~ . 
Dee 
de 

з 

У 

Mo 1 


та 
ita] 


= 
г.) 


Да + 





— | 
— 








Ее: 
ie d 





, : ДЕ 
. . bet zip be 
: ! rr . > s 
Е. 122) ed 
4 "ev, - 6 |- 
i -tte 44: . 
- wu P 

-.%4... о - 

ЕЕ . 28 |! 

LE D - „ . 

.. | . е 


лее 
.) 





li. 
еее Е 
—— — 
cortege 


t : 


кош кч 


*»b- 
D 
114 
ји 


l œj o = 
Pac 
zu 





A wA ——ÍÀ IO a A qu meer чий? cm p a m > em A AAA A ee LIP A md о 0 As. 





У | TECHNIQUE 
НИ: NE: ре ee o a L CONCENTRATION _ 


@----- q-- -— 4 · - .. -4 -: .--4 ZELE | ar --4 nn 4 —4-- = 4 И 


#75) 


| 


zer 


ves for 5-5 9.5mm polari 


y 
= 


са а а сте e 


[xi 


66 


р } ) Б | ----- )|-- 


—  —Á =... => ә 


„АМ ИМЕ : RESPONSF can" EXP 
2 3 








4 
: ‘| SUTWIDTH __ 





MINS 


- = еа» а Ф 
. 
• . 


5 | 








: Е : 
i : а =: . О. 
i - а а. ЗЕ ый ~ = ee a a - ree ИГ . E - 
i ! 4 4 ue к e . . 
Ие” | | | Е Е: 
— + - хе - = = еа =. уф n == e cm. «a -ea > тон е. = — а жиы 
| І NES X р, eet s 
ee а = ..4.. 2 2 ams в 15] 1“ “а .. езе... її = be г жә Ә----.. "ENS 
: : 2 Etc * ES de . «Un Coe = = Het трае: - em 4 dS eei 4 -... S .%-----. 4 
„| ар века ale Senn ni buceo ee Eu cce с TUE MEHR 
рр р О Дева re > ‚т 1р ари ӨЛШЕ irs Шү: 
е в ив . жа i. Lo.» есе и en + е ET ui - . . * à . ? LM 5 . ..ГИ . 5 E Ф р + 0 v 1. арии озш: A ани a | Ө кагу (At: 
500 Bi 2000. 1800 1000 600 


TECHNIQUE 


PATHLENGIY ©. Е А + 
DATE Ki An UNICA 


4 
CONCENTRATION REF No ла / А Ne-order No 


REFERENCE c — n — 1... Reerder No 


> = emm чин 
.- -Ф--.-.-Э---.. -.9——. B- - -H . -8— ———9 ———9—-- .9- —-— 39 ———- 9 —.–- 9 —–—-. 5—-–: –--9–--------В--—-- BD -- -.---- - - – 











è 
‘ 
‘ 
DAVE SS i 41: 
я NS . 2 » 
. ера роса е. сега 7° 
Фи пано в HEDS 
te. ne 
у е 
е 
~ 
. 


KRS-5 9.5mm polarize 


(#74) 


imi 


FOr 


erformance curves 


Б 


Figure 323 


67 








APPENDIX C 


R/V POINT SUR 


Description : 


R/V POINT SUR general specifications and dimensions. 


Conditions : 


Data extracted from R/V POINT SUR Cruise Planning Manual. 


Notes : 


None 
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APPENDIX D 
IMAGE FILES 


ЩЕ LIST OF THE SELECTED IMAGE FILES (*.PTW) 


Description : 


List of the *.PTW files showing their size, last modification date (creation date). 


Conditions : 


The list was generated directly from the Norton File Manager. 


Notes : 


None. 


ДВ 


Name 

<> 
NSA1001.PTW 
NSA1002.PTW 
NSA1004.PTW 
NSB1000.PTW 
NSB1002.PTW 
NSB1003.PTW 
NSB1005.PTW 
NSB1006.PTW 
NSB1008.PTW 
NSB1009.PTW 
NSB1011.PTW 
NSB1012.PTW 
NSB1015.PTW 
NSB1016.PTW 
NSB1018.PTW 
NSB1019.PTW 
NSCA0912.PTW 
NSCA0913.PTW 
NSCA0915.PTW 
NSCA0916.PTW 
NSCA0918.PTW 
NSCA0919. PTW 
NSCA0927.PTW 


ООО = 


NSCA0930.PTW 
МЕСАООВ РТ... 
NSCA0933.PTW 
NSCA0934.PTW 
NSOANUSSOET — == 
NSCA0937.PTW = 
NSCA0938.PTW 
NSCA0939.PTW 
NSCA0940.PTW - · 
NSCA0941.PTW -- 
NSCA0942.PTW 
NSCA0943.PTW 
NSCA0944.PTW 
NSCA0945.PTW -- 
NSCA0946.PTW 
NSCA0947.PTW 
NSCA0948.PTW 
NSCA0949.PTW 
NSCA0950.PTW 
NSCA0951.PTW 
NSCA0952.PTW 
NSCA0953.PTW 
NSCA0954.PTW 
NSCA0955.PTW 
NSCA0956.PTW 
NSCA0957.PTW — 
NSCA0958.PTW 
NSCA0959.PTW — 


NSCAO0902:PTW o e ee 
NSCA0963.PTW . mins 
NSCAQ0964.PTW `1 —- 


Size 


3,604,480 
1,703,936 
3,538,944 
1,703,936 
3,604,480 
1,703,936 
3,473,408 
1,769,472 
3,538,944 
1,835,008 
3,670,016 
1,703,936 
3,538,944 
1,769,472 
3,604,480 
1,703,936 
3,473,408 
1,703,936 
3,473,408 
1,703,936 
3,538,944 
1,835,008 
3,801,088 


1,835,008 : 


3922100 
1,769,472 
3,997,696 
1,900,544 


--3,670,016 


1,835,008 
3,473,408 
1,769,472 


3,473,408 - 


1,769,472 
3,407,872 
1,769,472 
3,538,944 
1,769,472 
3,473,408 
1,769,472 
3,604,480 
1,835,008 
3,932,160 
1,835,008 
3,407,872 
1,769,472 
3,473,408 
1,769,472 
3,473,408 
1,835,008 


. 3,604,480 " 


1,835,008 
3,538,944 
1,835,008 
3,604,480 


И 


Created On 
1/1/01 
3/9/98 
4/10/96 
4/10/96 
4/10/96 
4/10/96 
4/10/96 
4/10/96 
4/10/96 
4/10/96 
4/10/96 
3/9/98 
4/10/96 
4/10/96 
4/10/96 
4/10/96 
4/10/96 
4/9/96 
4/9/96 


- 4/9/96 


4/9/96 
4/9/96 
4/9/96 


4/9/96 
4/9/96 


4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 | 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 * 
4/9/96 
4/9/96 
4/9/96 


-201:50 PM ===> co EE 


11:22:30 AM 
8:27:24 AM 


... 8:28:48 AM 


8:32:54 АМ 

8:35:02 АМ 

8:35:00 АМ 

8:36:48 AM 

8:37:20 AM 

8:38:34 AM 

8:39:44 AM 

11:51:48 AM 

8:41:42 AM 

8:49:46 AM 

8:49:48 AM 

8:52:00 AM 

8:52:14 AM 

12:13:28 PM 

12:14:32 PM 

12:18:00 PM 

12:20:08 PM 

12:23:56 PM 

12:25:34 PM 

12:41:04 PM 

12:41:38 РМ — c P с сот 
12:44:30 PM 

12:44:48 PM 
12:46:30 PM 

12:47:06 PM 

12:51.06 PM -= ап 
12:51:20 PM 
1:03:50 PM 
1:05:08 PM 
1:07:06 PM ` 
1:08:42 PM _ 
1:10:34 PM 
1:11:38 PM 
1:15:00 PM 
1:15:40 РМ. 
1:19:44 РМ 
1:20:12 РМ 
1:27:42 РМ 
1:29:00 РМ 
1:30:28 РМ 
1:31:24 РМ . · 

1:3424 PM. -- ИМ . fc 
1:35:26 PM — | 

1:37:30 РМ 

1:38:38 РМ . | 
1:40:30 PM 2 zz NR 
1:4118РМ---:----- 2 сос 
132:56PM . . OE 
1:53:18 РМ -.:: _.. --- - E E - 
1:58:56 PM "m E 
1:59:54 РМ - Е А ISO) 


— m. 


= ama E 


Ф a m o w- 


Name 


NSCA0965.PTW 
NSCA0966.PTW 
NSCA0967.PTW 
NSCA0968.PTW 
NSCA0969.PTW 
NSCA0970.PTW 
NSCA0971.PTW 
NSCA0972.PTW 
NSCA0973.PTW 
NSCA0974.PTW 
NSCA0975.PTW 
NSCA0976.PTW 
NSCA0977.PTW 
NSCA0978.PTW 
NSCA0979.PTW 
NSCA0980.PTW 
NSCA0981.PTW 
NSCA0982.PTW 
NSCA0983.PTW 
NSCA0988.PTW 
NSCA0989.PTW 
NSCA0990.PTW 
NSCA0991.PTW | 
NSCA0993. TW 7c 7 20 
NSCA0994.PTW 
NSCA0996.PTW 
NSCA0997.PTW 
NSCA0999.PTW 
NSCA1086.PTW 


NSCA1087.PTW _ 


ee ee 


NSCA1089.PTW 
NSCA1090.PTW 
NSCA1092.PTW — 
NSCA1093.PTW 
NSCA1095.PTW 
NSCA1096.PTW 
NSCA1098.PTW 


NSCA1099.PTW__ 


NSCA900.PTW 
NSCA902.PTW 
NSCA903.PTW 
NSCA905.PTW 
NSCAS06.PTW 
NSCA911.PTW 
NSCA912.PTW 
NSCA914.PTW 
NSCA915.PTW 
NSCA917.PTW 
NSCA918.PTW 
NSCA920.PTW 
NSCA921.PTW 
NSCA923.PTW 


-. ба om ow Qi —À 


- come eee tee ee ome a 
- 


МАО ТИ Ре а аи 


NSCA926.PTW 


NSCA9Z7PTW — 
NSCA929.PTW = 7 


1 — = omm nn 


meee a ee ee ee 


--... 


of we es се 


Size Created On 
1,835,008 4/9/96 2:03:40 PM 
3,473,408 4/9/96 2:38:14 PM 
1,835,008 4/9/96 2:39:56 PM 
3,801,088 4/9/96 2:41:52 PM 
1,703,936 4/9/96 2:42:24 PM 
3,801,088 4/9/96 2:44:26 PM 
1,703,936 4/9/96 2:45:48 PM 
3,670,016 4/9/96 2:46:40 PM 
1,835,008 4/9/96 2:47:18 PM 
3,866,624 4/9/96 2:49:44 PM 
1,835,008 4/9/96 2:50:40 PM 
3,997,696 4/9/96 2:52:00 PM · 
1,900,544 4/9/96 2:52:28 PM 
3,735,552 4/9/96 3:04:40 PM 
1,769,472 4/9/96 3:05:24 PM 
3,735,552 4/9/96 3:06:30 PM 
1,703,936 4/9/96 3:07:14 PM 
3,538,944 4/9/96 3:09:34 PM 
1,835,008 4/9/96 3:10:14 PM 
3,670,016 4/9/96 3:18:04 PM 
1,769,472 4/9/96 3:19:38 PM 
3,670,016 4/9/96 3:44:46 PM 
1,835,008 4/9/96 3:45:02 PM 
3,735,552 4/9/96 8:50:22 5 т 
1,835,008 4/9/96 3:51:46 PM 
4,194,304 4/9/96 3:53:44 PM 
1,835,008 4/9/96 3:55:00 PM 
3,801,088 4/9/96 3:56:48 PM 
3,604,480 4/10/96  8&:1012АМ_______ _ ____ 
1,835,008 4/10/96 8:10:20 АМ 
3,932,160 4/10/96 8:17:14 АМ 
1,769,472 4/10/96 8:18:04 AM "ЖИЕ 
4,063,232 4/10/96 8:19:20 АМ = МЕ 
1,769,472 4/10/96 8:20:32 AM 
3,473,408 4/10/96 8:21:26 AM 
1,769,472 4/10/96 8:22:28 AM 
3,538,944 4/10/96 8:24:36 AM 
1,703,936 4/10/96 8:24:32 AM — Е дыз 
1,769,472 4/9/96 3:58:32 РМ 
3,670,016 4/9/96 4:01:06 РМ 
1,769,472 4/9/96 4:01:22 РМ 
3,604,480 4/9/96 4:03:14 РМ 
1,900,544 4/9/96 4:04:52 РМ 
3,538,944 4/9/96 4:20:16 РМ 
1,835,008 4/9/96 4:21:28 РМ | 
- 3,473,408 4/9/96 4:23:32 PM... 
1,835,008 4/9/96 4:23:32 PM 
3,670,016 4/9/96 4:26:08 PM 
1,769,472 4/9/96 4:26:08 PM — СТ 
3,604,480 4/9/96 428 12 Ро. xm 
1,769,472 4/9/96 ` 428:20PM ^ 7 
' 3,670,016 4/9/96 4:31:30 РМ. eo oen 
2 1,769,472 419196 -- : AOE PM бе 
3,735,552 49906 ____ 4ЗЗМАРМ ____________- 
1,835,008 4/9/96 4:340РМ..--2-..--- тв 
- 3,538,044 20 749196 7 4:51:06 PM REL CDI EEUU 


Name 

NSCA930.PTW 
NSCA932.PTW 
NSCA933.PTW 
NSCA935.PTW 
NSCA936.PTW 
NSCA938.PTW 
NSCA939.PTW 
NSCA941.PTW 
NSCA942.PTW 
NSCA944.PTW 
NSCA945.PTW 
NSCA947.PTW 
NSCA948.PTW 
NSCA958.PTW 
NSCA959. PTW 
NSCA961.PTW 
NSCA962.PTW 
NSCA964.PTW 
NSCA96S.PTW 
NSCA967.PTW 
NSCA968.PTW 
NSCA971.PTW 
NSCA972.PTW 
NSCA974.PTW 
NSCA97S.PTW 
NSCA977.PTW 
NSCA978.PTW 
NSCA980.PTW 


NSCA981.PTW.. . 


NSCA983.PTW 


Size 

1,769,472 
3,735,552 
1,900,544 
3,604,480 
1,769,472 
3,538,944 
1,703,936 
3,670,016 
1,835,008 
3,670,016 
1,769,472 
3,735,532 
1,769,472 
3,735,552 
1,769,472 
3,670,016 
1,703,936 
3,801,088 
3,407,872 
3,473,408 
1,703,936 
3,538,044 
1,835,008 


3,604,480 -- 


1,835,008 
3,670,016 
1,703,936 
3,670,016 


- 1,769,472. -.- 


3,538,944 


Created On 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 
4/9/96 


4/9/96 


4/9/96 
4/9/96 
4/9/96 - 
4/9/96 
4/9/96 
4/9/96 
4/9/96 


- 4/9/96 


4/9/96 


4:51:20 PM 
4:54:02 PM 
4:54:12 PM 
4:55:38 PM 
4:56:36 PM 
4:57:44 PM 
4:58:46 PM 
5:00:16 PM 
5:01:22 PM 
5:13:34 PM 
5:15:00 PM 
5:16:28 PM 


5:16:24 PM 


5:30:28 PM 
5:31:20 PM 
5:33:16 PM 
5:34:34 PM 
5:43:36 PM 
5:43:36 PM 
5:45:16 PM 
5:46:40 PM 
5:48:08 PM 
5:49:40 PM 


ОР. 


5:51:20 РМ 
5:53:56 РМ 
3:53:52 PM 
3:55:56 PM 


5:56:02 PM -—-- 8 


5:57:12 PM 


2. LIST OF THE BASIC IMAGE FILES (*.PTE) 


Description : 


List of the *.PTE files showing their size, last modification date. 


Conditions : 


The list was generated directly from the Norton File Manager. 


Notes : 


None. 


7? 


Name 

BASEOI.PTE 
BASE02.PTE 
BASEO3.PTE 
BASE04.PTE 
BASEOS.PTE 
BASE06.PTE 
ВА5Е07.РТЕ 
BASE08.PTE 
BASEO9.PTE 
BASEIO.PTE 
ВАЗЕПРТЕ 
BASEI2.PTE 
BASEI3.PTE 
BASE14.PTE 
BASEIS.PTE 
BASEI6.PTE 
ВА5Е17.РТЕ 
BASEIS.PTE 
BASEI9.PTE 
BASE20.PTE 
BASE21 PTE 
BASE22. PTE 
BASE23.PTE 
BASE24.PTE 
BASE25. PTE 
BASE26.PTE 
BASE27.PTE 
BASE23.PTE 


BASE29PTE ” zer: 


ВАЗЕЗОРТЕ 2 
BASE31.PTE 
BASE32.PTE 
BASE33.PTE 
BASE34.PTE 
BASE35.PTE 
BASE36.PTE 
BASE37.PTE 
BASE38.PTE 


BASE40.PTE 
BASEALPTE - - 
BASE42 PTE 
ВА5Е4ЗРТЕ 
BASE44.PTE 
BASE45 PTE 
BASE46 PTE 


ВАЗЕД7РТЕ---------- 7” 


BASE48 PTE 
BASE49 PTE 
BASESO.PTE 


PTE File 
PTE File 
PTE File 
PTEFile 
PTE File 
PTE File 
PTE File 
PTE File 
PTE File 
PTE File 
PTE File 
PTE File 
PTEFile 
PTE File 
PTE File 
PTEFile 
PTE File 
PTE File 
PTE File 
PTE File 
PTE File 


PTEFile . 


PTE File 
PTEFile 
PTE File 
PTEFile 
PTE File 
PTEFile 


------- | 


PTE File 


PTE File —- — 


PTE File 
PTEFile 
PTEFile 
PTE File 
PTEFile 


Pie File ли 


PTE File 
PTE File 


PTE Fie... 


PTE File 


Pre File = 


PTE File 
PTE File 
PTE File 
PTE File 
PTE File 


РТЕЕе ^ | 


PTE File 
PTE File 
PTE File 
PTE File 
PTE File 
PTE File 
PTE File _ 
PTE File 
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Size 

4,896254 
4,896,254 
5,025,070 
5,347,110 
5,347,110 
5,540,334 
5,089,478 
4,831,846 
4,831,846 
4,767,438 
4,896,254 
4,896254 
5,089,478 
5,411,518 
4,831,846 
4,960,662 
4,896,254 
5,025,070 
5,025,070 
5,089,478 
4,896,254 
5,089,478 
5,089,478 
5,089,478 
5.282,702 
3,475,926 
5,218,294 
5,025,070 


4,960,662 — 


5,025,070 
5,089,478 
5,153,836 
5,540,334 
5,153,886 
5,025,070 


5,089,478. 


5,025,070 
4,896,254 


5,025,070 ЕЕ 


5,025,070 
5,025,070 
5,153,886 
4,960,662 
5,282,702 
5,089,478 
4,831,846 
5,089,478 
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5,153,886 
5,089,478 
5,025,070 
5,797,966 
4,767,438 
4,960,662 
5,025,070 
4,960,662 
5,089,478 
5,025,070 
5,282,702 
5,411,518 
4,896,254 
4,896,254 
5,025,070 


4,960,662 .. 
" 5,025,070 -- ---- 
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5,089,478 
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5/24/98 6:08:38 PM 

5/24/98 6:24:48 PM 
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E 7:13:36PM — 1 _ 
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APPENDIX E 


METEOROLOGICAL DATA 


|. METOC FILES 


Description : 


This Appendix presents the description of METOCI and METOC2 files. Their recorded 


parameters are listed and a sample taken directly from the METOC files is shown. 


Conditions : 


Information provided by NPS-Meteorological Group. 


Notes : 


None. 
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STORAGE FORMAT FOR METOCO01 (DASC) DATA 


The data is stored in a set of files with the names "eopace.*" and "30sec.*" The extension 
for the file names is the Julian date the file was started on. A new pair of files is started every day 
at 00:00 GMT. The time data in the files are all GMT. The eopace files have the data averaged 
over 10 minutes. The 30sec files have the GPS data averaged over 30 seconds and the 
last set of data from the Campbell. In the 30sec file, the Campbell data are not averaged. The 
data format for both files is the same. Each file has a header which tells when the file was started. 
Each file also has a footer that tells when the file was finished. Some files may also have 
additional start and stop times interspersed with the data, depending on whether data logging was 
turned off or the program was stopped before the end of the day. The data are arranged in fields 
separated by spaces across a line of text. A carriage return character and a line feed character are 
at the end of each line. Each new reading is on a new line. The fields are in the order below. The 
numbers before the fields are not in the data. They serve to make identifying the fields easier. 

01 number of readings in this average, nn 
02 year, yyyy 
03 Julian date, ddd 
1—— — поште бе minutes; hhmm- — — ——— —————————————————— Ы 
05 seconds, ss 
06 Relative Wind Speed, m/s, mmm 
07 Relative Wind Direction, degrees, ddd 
08 Ship Speed, knots, at present always 0 
09 Ship Direction, degrees, at present always 090 
10 True Wind Speed, m/s, mmm 
11 True Wind Direction, degrees, ddd 
12 T air, degrees C, tt.t 
13 RH, per cent, pp 
14 Pressure, millibars, pppp.p | 
15 Sea Surface Temperature, degrees C, tt.t 
16 GPS time, hhmmss 
17 GPS Latitude, ddmm.m 
18 GPS North/South Indicator, either 'N' or 'S' 
19 GPS Longitude, dddmm.m 
20 GPS East/West Indicator, either 'E' or 'W' 
21 GPS Speed Over Ground, knots, kk.k 
22 GPS Course Over Ground, degrees, ddd.d 
23 GPS Antenna Height, meters, mm.m 


Here is a sample of the data: 
30 1996 99 1103 28 5 009 0 090 3 323 123 91 1013.9 
20.3 110838 3306.2 N 11806.0 W 0.0 90.0 189.6 
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STORAGE FORMAT FOR METOC02 DATA 


The METOC02 system had two additional temperature sensors and two additional RH sensors. The 
data are arranged slightly differently than the METOCO1 system. The data is stored in a set of files with the 
names "eopace.*" and "30sec.*". The extension for the file names is the Julian date the file was started on. A 
new pair of files is started every day at 00:00 GMT. The time data in the files are all GMT. The eopace files 
have the data averaged over 10 minutes. 

The 30sec files have the GPS data averaged over 30 seconds and the last set of data from the Campbell. In 
the 30sec file, the Campbell data are not averaged. The data format for both files is the same. Each file has a 
header which tells when the file was started. Each file also has a footer that tells when the file was finished. 
Some files may also have additional start and stop times interspersed with the data, depending on whether 
data logging was turned off or the program was stopped before the end of the day.The data are arranged in 
fields separated by spaces across a line of text. A carriage return character and a line feed character are at the 
end of each line. Each new reading is on a new line. The fields are in the order below. The numbers before 
the fields are not in the data. They serve to make identifying the fields easier. 

01 number of readings in this average, nn 

02 year, yyyy 

03 Julian date, ddd 

04 hours & minutes, hhmm 

05 seconds, ss 

— — — — -706Relative Wind Speed, m/s,;inmn | aS 

07 Relative Wind Direction, degrees, ddd 

08 Ship Speed, knots, at present always 0 

09 Ship Direction, degrees, at present always 090 

10 True Wind Speed, m/s, mmm 

11 True Wind Direction, degrees, ddd 

12 T air, degrees C, tt.t 

13 КН, регсеп рр (0 

14 Pressure, millibars, рррр.р 

15 Sea Surface Temperature, degrees C, tt.t 

16 T air A, degrees C, tt.t 

17 RH A, per cent, pp 

18 T air B, degrees C, tt.t 

19 RH B, per cent, pp 

20 delimiter between Campbell data and GPS data, "***" 

21 GPS time, hhmmss 

22 GPS Latitude, ddmm.m 

23 GPS North/South Indicator, either 'N' or 'S' 

24 GPS Longitude, dddmm.m 

25 GPS East/West Indicator, either 'E’ or "А 

26 GPS Speed Over Ground, knots, kk.k 

27 GPS Course Over Ground, degrees, ddd.d 

28 GPS Antenna Height, meters, mm.m 


Here is a sample of the data : 
18 1996 99 0319 28 2 090 0 090 3 263 122 96 1012.3 16.4 
12.2 96.9 12.7 95.7*** 031908 3322.3 N11836.9 W 5.0 1263 25.7 


METOC 1 FORMAT 


Data logging started at Tue Apr 09 00:00:00 1996 


30 


21.2 000024 3311.2 


30 


21.2 000054 3311.2 
1996 
21.2 000124 3311.1 


30 


30 


21.2 000154 3311.1 


30 


21.2 000224 3311.1 


30 


21.2 000254 3311.1 


30 


21.2 000324 3311.0 


30 


21.2 000354 3311.0 


30 


21.2 000424 3311.0 


1996 


1996 


1996 


1996 


1996 


1996 


1996 


1996 


99 2355 44 2 156 O 090 
N 11728.2 М 63 1243 56.8 
2355 44 2 156 0 090 


INZEIN2 5:29 9996: 3.123. 9 525 
2355 44 2 156 0 090 
М 11728.1 \ 63 1235 42.0 


2355 44 2156 ШО 090 
М 117280 ХУ 63 123.0 26.9 
2355 44 2 жшс 0 090 


N 117280 W 6.4 121.6 99 


Data logging started at Tue Apr 09 00:00:00 1996 


30 
16.8 
30 
16.8 
30 
16.8 
Zo 
16.8 
30 
16.6 
30 
16.6 
30 
16.6 
30 
16.6 
30 
16.6 


1996 
0.0 
1996 
0.0 
1996 
0.0 
1996 
0.0 
1996 
0.0 
1996 
0.0 
1996 
0.0 
1996 
0.0 
1996 
0.0 


0.0 0.0 0.0*** 000035 
09 2352 32 2 
0.0 0.0 0.0*** 000105 
20235220057 2 
00 00 0.0*** 000135 
МЫР 233250532 2 
0.0 0.0 0.0*** 000205 
100 0002 8 2 

0.00 0.0*** 000235 
100 0002 8 2 

0.0 0.0*** 000305 
100 0002 8 2 


0.0 


0.0 


2355 44 2 die, 0 090 

N 11727.9 W 63 1205 -42 
2355 44 2 156 O 090 

N 11727.9 W 5.9 129.6 -10.7 
2355 44 2 156 0 090 

N 117279 W 4.1 185.6 -6.7 
2355 44 2 156 0 090 

N 11727.9 W 4.1 233.5 9.0 
METOC 2 FORMAT 
Ао 2 138 0 090 


138 0 090 


138 0 050 
138 0 090 
149 O 090 
149 0 090 


149 0 090 


3311.2 N 117282 


3311.1 М 11728.1 


3311.1 М 11728.1 


3311.1 М 11728.0 


3311.1 М 11728.0 


331140 TINI UP DES 


а лала X uA 


W 
5 


0.0 0.0 0.0*** 000335 3311.0 N 11727.9 W 
100 0002 8 2 149 0 090 5 
00 00 0.0*** 000405 3311.0 М 11727.9 W 


100 0002 8 2 
0.0 00 0.0%% 000435 


149 0 090 


62 


D 


заети 


294 


294 


294 


294 


294 


294 


294 


294 


294 


16.5 
124.3 
16.5 
124.2 
16.5 
| e 
16.5 
122.3 
16.2 
121.4 
16.2 
120.7 
16.2 
159.8 
16.2 
216.5 
16.2 
259.1 


81 1016.2 


81 1016.2 


81 1016.2 


81 1016.2 


81 1016.2 


81 1016.2 


81 1016.2 


81 1016.2 


81 1016.2 


75 1016.9 
46.3 
75 1016.9 
40.7 
75 1016.9 
29.4 
75 1016.9 
14.0 
76 1016.9 
-2.1 
76 1016.9 
-13.8 
76 1016.9 
-17.5 
76 1016.9 
-10.2 
ооо. 
8.3 


D METEOROLOGICAL CONDITIONS 


Description : 


This Appendix presents the experimental and the calculated conditions for all 70 analyzed 
files. The data 15 presented in 2 tables where each column represents a parameter and each line 
represents a file. File number, date, recording time, GPS data, angular positioning and 
atmospheric transmittance were added to the METOC data. 


Conditions : 


Meteorological data presented in this Appendix represent the average value from 
METOC1 and METOC2 during the record time interval of each file (image sequence). Range and 
angular information were obtained by calculation using the ship GPS and sensor location 
information. Atmospheric transmittance was calculated using MODTRAN under conditions 
described in Chapter III item A.4.c. 


Notes : 


Day column shows 9 or 10 corresponding to April 9 and April 10, 1996. 

Values in "time" column are elapsed time given in seconds starting at 00h:00min (local) of 
the corresponding day. 

Latitude (Lat) and longitude (Lon) are presented in minutes and the ship GPS antenna 
height (H) is presented in meters. 

Heading and Bearing information are referred to the True North. 

The range and angular information (Range, Theta, Phi) are spherical coordinates referred 
to the camera coordinate system: 


- Origin : Camera position 

- x axis : Oriented to true North 

- y axis : Oriented to true West = 
- Z axis : Up 


True North 
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APPENDIX F 


SHIP SKIN TEMPERATURE 


Description : 


This Appendix presents the position of the thermistors used to measure the ship skin 
temperature during the experiment. It also presents the format used to store the skin temperature 
data. 


Conditions : 
The data was collected by the NPS-NACIT group using 16 thermistors in different 
positions on the ship surfaces. 
Notes : 


The first column of the skin temperature data is the elapsed time in seconds from the initial 
recording time. This initial time is recorded at the top of the file. 
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SKIN TEMPERATURE DATA FORMAT 


april9. 1 
Mon,Day,Year- 4 9 1996 
Hour,Min,Sec= 20 8 27 
volts,istart,istop,nsam,iref 1.49939 1 16 100 1 

763707 25.001 14.458 17.209 13.710 14.943 16.570 13.946 14.468 25.001 15.562 14.130 
31.901 15.310 27.986 15.783 15.153 

763727 25.001 14.451 17.108 13.704 14.938 16.585 13.948 14.457 25.001 15.561 14.115 
585826 15.315 27.927 15.780 15.149 

763747 25.001 14.448 17.189 13.692 14.928 16.573 13.952 14.463 25.001 15.553 14.110 
31.801 15.299 27.878 15.776 15.148 

763767 25.001 14.435 17.225 13.679 14.922 16.579 13.944 14.450 25.001 15.547 14.099 
31.789 15.295 27.852 15.771 15.141 

763787 25.001 14.424 17.226 13.666 14.916 16.584 13.927 14.436 25.001 15.539 14.089 
31.769 15.311 27.839 15.769 15.148 

763807 25.001 14.422 17.225 13.660 14.911 16.591 13.922 14.425 25.001 15.546 14.084 
31.744 15.308 27.816 15.769 15.134 

763827 25.001 14.420 17.226 13.649 14.908 16.599 13.922 14.411 25.001 15.537 14.071 
31.701 15.285 27.814 15.766 15.135 

763847 25.001 14.420 17.221 13.650 14.907 16.608 13.922 14.405 25.001 15.550 14.066 
31.628 15.266 27.744 15.764 15.129 


763867 = ELAPSED TIME IN SECONDS FROM INITIAL TIME 


25.001 14.417 17.224 13.648 14.905 16.610 13.920 14.412 25.001 15.562 14.070 31.528 
15.233 27.655 15.757 15.129 = TEMPERATURE FROM EACH THERMISTOR IN C 
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APPENDIX G 


AGA 780 LABORATORY CALIBRATION CURVES 


Description : 
This Appendix presents the laboratory calibration of the camera. 


Conditions : 


The calibration was done in MAR/21/96 by E.C. Crittenden for the NOSC LW Head with 
internal polarizers. 


Notes : 


Constant conversion for use in CEDIP or CATS software is provided. 

Refer to the AGA780 manual [Ref. 8] for further information. 

The calibrations is done for each filter setting by recording the thermal level measured of a 
blackbody for a set of different temperatures. 

The calibration curve is then obtained by a curve fitting process. The equation used for the 
curve fitting has a set of constants. Depending on the acquisition software (AGEMA-CATS or 
CEDIP PTRWIN), a different set of constants is determined by the curve fitting process. In this 
Appendix are presented both sets of constants as well as the conversion equations. 
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Calibration of the NOSC LW Head with internal polarizers. 
Ван кеси Е С. Crittenden. 

Using CATS calibration constants. 

Ш ЕЕ ос data are in "Split Field" notebook (green 
hardback) 

Copies are included with this. 

CATS constants are A, B,OS and C, with C always = 1. 
CEDIP constants are R,B, V, e, and F with F always = 1. 


The constants obtained are; 


Pos A B OS 
4 1675 1348 39 
5 1675 1348 39 
6 6808 1511 12 


The corresponding values for the CEDIP constants are obtained 


re ee ee: ш — + nn nn —— 9 ооз ч ни — - ---- -- --- = ——— amacaa 


from the equations on the following page. 


пр APPS — PPP et ee oe — = 1.. sæ: omo mom — ---- .--.. - --. -- - - т... әм. nn nn -- -- — ———— = .—_ ——— –— - 
- =- -- 
— = - = --— — = ——— eee — s 9 À AA — — M o  ———MAÀÀÀ o — шл rn - 
= = — = - m = -= et AA === — mI 
.. --. -- 
` —— | - — — ——— —— ———————P———ÁÓÓ (Po m 
.. 
— — ээ» re itae ви = - -- - - е -- ------- --- > rm o —À = +o o --- 
== а == — — M өлке 
a m -- — 5 а => =- u rn ee ee nn an ст ъв. + 
-. = ---- - = == - - >. 
= ~ = - x 
+ — — - >. y cer 
= === we === - — 
qe A ee —— nen PP то -— - - ——— — - — > ame 
> - — -- — 
— . .. 
—— o ee с ыш лалы... 2 ананна —————— 
. а а» п 4 mum — u —————————————————————Ó—Á——————— 
— >. == .- - -- 
- - — — — ~ == —ÁÀ ~ = — ———— BR - 
. 
> - 


7 
О 
Розшоп 4 


x 


) Head. Lw 


, 3/21/96 


і 
i 
i 


¡osition 6 


|. nr ng 


| 


(N 


ы [заа 
| 


© 


Internal Pol 
u 
Bc 
И 
60 


OSO 
arizers 
Li 
eee ep | 
Tl 
INE 
- 40 
—70 
CAVITY BLACK! BODY TEMPERATURE 


| 
| 


NRAD 
өз: 
—80 
—90 
| 
Position 5 


c Mae 
Ен 

- 110 
+ 


ES O E 
O 
N 
J = 
вр | 
a 


> бо. = © D = 5 < са 


2 (3 гомаузи BAT IVWSIHL- от 


.— een = ---.- -- --.. == am eoo omar o -- 


——-.— ern m = е mann 
нее o mu ав = — 
В . . > - O a a 
—-——— ць 

—— — — un «т — —————— MÀ — M — а ош ee oe — 
. - - AA ee ee eee = —— 


| a 
=... =- .- — - - © = o. —. 
7 = une el em amem от - - ------ - 
FR eu une um = ee 5 Е --- -—— oo 949 9m жо» ... а cc C аа m ne 
a 20 € —Ó и as <oe ees oe —— — ———— oar ои s sem 
у — - - - — ~ -- cx 
TO у ро Sw eS ње и срне P = А 
+ . - —À— as Sat He 9e, P 
E м - -- — e — a æ -<= -- 
— —= = nn. - - .- - > ey ee T қ, Dr ER NE EDAD RO --.-...... 
ые...” ee тез : + . . = . : iras E А, S api E ^ RE Tm A [ыы m rm = -— = 
| SA - o s= ...-- 
ste Ды = А --. 1 ma s asime ao apet 
FA AX a 
æ . ... -- - Maur uu o D Ре rre . =. ас = 
-. - - -. 5 — —À e. - а 
.- EM — —MUnáóáÀ € "— тр А - E - - - > o Пати СА 4 
E EIS маа a .--. --- -а - .- - I A Е а. аана —  — > - ow > Е - и 
- a я . .. = ЕЕ --...- «важна Ее. - =.. ----- MÓÓ— M ——— e e 
. .. Wa) (m. ми а о а zi xi B "ees TT TM — À — —À M « 
Wee ee PE SSS Ири Мы — — sty em a ee en en m ee . г 
— — — — - г - 
— --- — — — x A m -- -— — -— — - — — - = Е 
— - — 
ee ee eee = — — — — 
= -- - — — om — — = 
— - - — ] - 
-- петљи тв њи = — te te om — --> 
. ~. е 
— = a - =. = = - 


o — PQ 0 = 
dL re eee 


у E: CDP-CNST | CT 5/16/96 


"SUMMARY OF VALUES TO ENTER IN CEDIP: 
Numericel velues below ere exemples: 














AGEMA Thermal level (TL): (C= 1 изиецу] 
A= 3708 From AGEMA celibretion curve fit | 
А 
B= 1447 СЕОТР В = АбЕМА В ІШЕ - % (05) 
C*EXP(E/T) - 4 
F= 1 Almost elweys 
v= 293 Called "Background Temp". This can be CEDIP Thermal level (UI): (Е = 1 usually) 
erbitrerily chosen, but will chenge s 
end hence R. Set V to 273 + the temp R R 
in C. most important in the imege. (ur - в ---- + [1]? 
ехр(В/Т) ~ Е ехр(В/У) - Е 
0S = 10.4 0S from AGEMA calibration 
A 
еш 0.80 Calc, from formule at right—> e= ТОЛ ан = 0.7984 
(Two decimal pleces can be entered in CEDIP) OS*[exp(B/V) -F]*A 
R= 7149 Н = А/е 
| Ang A 
05 - 10.4 Celc. es e check, with formule et right —> 05 - 5 = 10,40 


е exp(B/V) - F 
ENTERING CALISRATION CONSTANTS: 


— = 


m" a г 8, end F. ere. веб in СЕОТР. under. > "Measure, Calibration, . New, (or Open)" и 
Torote O TT O E should normelly be Left et unity. 

e is set in CEDIP under "Measure,Emissivity" 

V is set in CEDIP, under "Meesure, Beckground" 


VERIFYING AND CORRECTING THE OFFSET: 


To verify the offset, meesure the temperature of a bleck body et the tempereture chosen for V. 
If en incorrect T is obteined, increesing V increases the offset (slides the whole curve uniformly upward) 
increasing M will lower the observed tempereture, 
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VERIFYING AND CORRECTING THE CALIBRATION CURVE SLOPE: ^ ^ 


This correction should be rarely needed. : 
Verify @ tempereture above, or below, the temperature chosen for V. : 
If that T is incorrect, incressing e increeses the célibration curve slope, rotating ebout the T = V value, 
Increasing e elso Lowers the offset, corresponding to rotetion of the curve about T = V, This offset 

lowers the curve to compensete for the new slope, so the curve still pesses through the point st T = V. 
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APPENDIX H 
ANALYSIS PROGRAMS 


1. LIST OF THE ANALYSIS PROGRAMS 


Description : 


This Appendix presents a table describing the type and the objective of all analysıs 
programs written in IDL. 


Conditions : 


None. 


Notes : 


None. 
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TABLE Hl 
ANALYSIS PROGRAMS 


elements by column, line and general 
Find the edges of the analysıs “box” for each frame 
Convert *.PTW files to *.PTE files 
О | 


Elements Аихшагу Generate “mask matrices” (1 for the object pixels and 0 otherwise) for : 
ship, sea, sky, ship above horizon and ship below horızon. 


Inclua Include data into the result matnx 
Leia_pte Read variables from basic files 


Linespace Create a vector of n numbers linearly spaced between two given values 





Meteoro Interpolate meteorological data to each frame time 

Modele Show a scaled ship model as surfaces or wireframe. Shows a ship model 
demonstration . Optionally the model can be drawn over a given image. 

Norma Normalize the image 

Pontes Run analysis programs in correct sequence 
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Mn parto | ore ener — 


1.This program is based on “draw_box.pro” program created by Mr. Jerry Lentz - Physics 
Department / 1994. 
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2 HOTSPOT TRACKING 


Description : 


This Appendix presents figures showing the hotspot program tracking feature. 


Conditions : 


Sample extracted from BASE18.PTE unpolarized case. The R/V POINT SUR in this file 
Is at range 5,134m, motionless near to the horizon line. However, in order to show the hotspot 
tracking capability in a cluttered scene (3 targets in this case), a small boat was chosen to be 
tracked on its hotspot. Therefore, in this Appendix, a sequence of 8 images is presented showing 
the tracking of a small boat moving from left to right at the bottom of the images. 


Notes : 
The hotspot position is indicated by a small square superimposed to the target image. The 


image sequence shows the hotspot tracking capability of the program HOTSPOT.PRO and 
possibility of using a cluttered image for analysis without losing the hotspot correctness. 
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HOTSPOT TRACKING IN A CLUTTERED SCENE (BASE18.PTE) 
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3: BOX POSITIONING 


Description : 


This Appendix presents figures showing the analysis box positioning. 


Conditions : 


Sample extracted from BASE18.PTE unpolarized case. The R/V POINT SUR in this file 
is at range 5,134m, motionless near to the horizon line. However, in order to show the box 
positioning capability in a cluttered scene (3 targets in this case), a small moving boat was chosen 
to be the target. Therefore, in this Appendix, a sequence of 8 images is presented showing the box 
positioning around a small boat moving from left to right at the bottom of the images. 


Notes: 


The analysis box is indicated by a square around the target image. The box positioning 
uses the hotspot position as its reference. 
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4. 


IMAGE SAMPLE 


Description : 


This Appendix presents one of the 5,508 analysed images. This is the frame #0 from the 


unpolarized case in BASEOI.PTE file (file 701). 


Conditions : 


From tables in Appendix E.2: 


File #01 : 


Notes: 


Day  : April 9 1996 

Time : 44,018 seconds (elapsed time from 00h:00min local time = 12h13min38sec) 
Latitude : 1959.20 minutes 
Longitude : 7034.80 minutes 

Ship GPS Antenna Height : -27.01 m 
Air Temperature : 15.30 C 

Relative Humidity : 72 % 

Air Pressure : 1018.25 mb 

Wind Speed : 6.5 m/s 

Sea Surface Temperature : 18.55 C 
Range : 998.52 m 

Theta : 1.62 rad 

Phi : 2.41 rad 

Heading : 99 degrees 

Bearing : 221.85 degrees 
Atmospheric Transmittance : 0.82 


Conditions defined in Appendix E.2. 
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5. 


COMBINED IMAGE / SURFACE / CONTOUR PLOT 


Description : 


This Appendix presents one ofthe 5,508 analysed images. This ıs the frame #0 from the 


unpolarized case in BASEO1.PTE file (file #01). The figure is in "show3" IDL format. This 
format shows a combined picture with image, surface and contour plots. 


Conditions : 


From tables in Appendix Е.2: 


File #01 : 


Notes: 


Day : April 9 1996 

Time : 44,018 seconds (elapsed time from 00h:00min local time = 12h13min38sec) 
Latitude : 1959.20 minutes 
Longitude : 7034.80 minutes 

Ship GPS Antenna Height : -27.01 m 
Air Temperature : 15.30 C 

Relatıve Humidity : 72 % 

Air Pressure : 1018.25 mb 

Wind Speed : 6.5 m/s 

Sea Surface Temperature : 18.55 C 
Range : 998.52 m 

Theta : 1.62 rad 

Phi : 2.41 rad 

Heading : 99 degrees 

Bearing : 221.85 degrees 
Atmospheric Transmittance : 0.82 


Conditions defined in Appendix E.2. 
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BASEOI.PTE / UNPOLARIZED CASE / FRAME #0 
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6. 


SURFACE PLOT 


Description : 


This Appendix presents one of the 5,508 analysed images. This is the frame #0 from the 


unpolarized case in BASEO1.PTE file (file #01). The image is shown in surface plot. 


Conditions : 


From tables in Appendix E.2: 


File #01 : 


Notes: 


Day : April 9 1996 

Time : 44,018 seconds (elapsed time from 00h:00min local time = 12h13min38sec) 
Latitude : 1959.20 minutes 
Longitude : 7034.80 minutes 

Ship GPS Antenna Height : -27.01 m 
Air Temperature : 15.30 C 

Relative Humidity : 72 % 

Air Pressure : 1018.25 mb 

Wind Speed : 6.5 m/s 

Sea Surface Temperature : 18.55 C 
Range : 998.52 m 

Theta : 1.62 rad 

Phi : 2.41 rad 

Heading : 99 degrees 

Bearing : 221.85 degrees 
Atmospheric Transmittance : 0.82 


Conditions defined in Appendix E.2. 
Main image elements are indicated in the picture. Observe the characteristic "wall shape' 


at the horizon line. 
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7: 


NORMALIZATION : PHASE 1 


Description : 


This Appendix presents one of the 5,508 analysed images during the normalization process 


used by the THRESHOLD.PRO program. This is the frame #0 from the unpolarized case in 
BASEO1.PTE file (file 201). The image is shown in surface plot. 


Conditions : 


From tables in Appendix E.2: 


File 201 : 


Notes : 


Day  : April 9 1996 

Time : 44,018 seconds (elapsed time from 00h:00min local time = 12h13min38sec) 
Latitude : 1959.20 minutes 
Longitude : 7034.80 minutes 

Ship GPS Antenna Height : -27.01 m 
Air Temperature : 15.30 C 

Relative Humidity : 72 % 

Air Pressure : 1018.25 mb 

Wind Speed : 6.5 m/s 

Sea Surface Temperature : 18.55 C 
Range : 998.52 m 

Theta : 1.62 rad 

Phi : 2.41 rad 

Heading : 99 degrees 

Bearing : 221.85 degrees 
Atmospheric Transmittance : 0.82 


Conditions defined in Appendix E.2. 
The main elements of the image are indicated. Compare with the surface plot of the same 


image before the normalization (Appendix H.6). After normalization (phase 1) the characteristic 
"elevation" seen in the sky thermal level surface plot disappeared. The depression caused by the 
target presence is visible. "Noise" level at the sky and sea image components are still significant. 


114 


ЕНСІН 


| HSVHd : NOLLVZTVIVIN?ION /04 HINV?IH / HSVO CIHZISIV'IOdN( / Ad TOUS Va 





ШЕ 


8. NORMALIZATION : PHASE 2 


Description : 


This Appendix presents one of the 5,508 analysed images during the normalization process 
(phase 2) used by the THRESHOLD.PRO program. This is the frame #0 from the unpolarized 
case in BASEOI.PTE file (file 201). The image is shown in surface plot. 


Conditions : 
From tables in Appendix E.2: 
File #01 : 


Day  : April 9 1996 

Time : 44,018 seconds (elapsed time from 00h:00min local time = 12h13min38sec) 
Latitude : 1959.20 minutes 
Longitude : 7034.80 minutes 

Ship GPS Antenna Height : -27.01 m 
Air Temperature : 15.30 C 

Relative Humidity : 72 96 

Air Pressure : 1018.25 mb 

Wind Speed : 6.5 m/s 

Sea Surface Temperature : 18.55 C 
Range : 998.52 m 

Theta : 1.62 rad 

Phi : 2.41 rad 

Heading : 99 degrees 

Bearing : 221.85 degrees 
Atmospheric Transmittance : 0.82 


Notes : 


Conditions defined in Appendix E.2. 

The main elements of the image are indicated. Compare with the surface plot of the 
normalization phase 1 (Appendix H.7). After normalization (phase 2) the "noise" level at the sky 
and sea image components are reduced and the background surface is basically flat. The 
depression caused by the target presence is still visible. The ship is now easily identified and there 
is no alteration in his shape. 
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9. 


IMAGE ELEMENTS 


Description : 


This Appendix presents one of the 5,508 analysed images during the image elements 


identification process done by the ELEMENTS.PRO program. This is the frame £0 from the 
unpolarized case in BASESO.PTE file (file 250). The image is shown in its full version with the 
hotspot and the analysis box indicated. Below the full image, the analysis box and the main 
element masks (ship, sea and sky) are shown. 


Conditions : 


From tables in Appendix E.2: 


File £50 : 


Notes : 


Day  : April 9 1996 

Time  : 62,187.5 seconds (elapsed time from 00h:00min local time — 17h31min8sec) 
Latitude : 1959.40 minutes 
Longitude : 7034.40 minutes 

Ship GPS Antenna Height : -19.44 m 
Aır Temperature : 15.30 C 

Relative Humidity : 69 96 

Air Pressure : 1016.40 mb 

Wind Speed : 4.0 m/s 

Sea Surface Temperature : 16.60 C 
Range : 373.21 m 

Theta : 1.71 rad 

Phi : 3.02 rad 

Heading : 68 degrees 

Bearing : 186.9 degrees 
Atmospheric Transmittance : 0.92 


The white color in the mask figures (ship, sea and sky) represents pixels with values of 1. 


The black color means pixel value of zero. 
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10. SHIP MODEL 


Description : 

This Appendix presents the ship model viewed from different angles. 
Conditions : 

Elevation angle for all picture is 2 degrees. 


Starting from the upper left picture, the azimuth angles are: 45, 0, 315, 135, 180 and 225 
degrees. 


Notes : 


None. 
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SHIP MODEL : WIREFRAME MODE 


Description : 


corresp 


This Appendix presents the ship model in wireframe mode superimposed to analysis box 
onding to the frame £0 of the unpolarized case of the BASESO.PTE file. 


Conditions : 


Model elevation and azimuth angle are 2 and 70 degrees respectively. 
For the Image, from tables in Appendix E.2: 


File 450 : 


Notes : 


Day :April9 1996 

Time :62,187.5 seconds (elapsed time from 00h:00min local time = 17h31min8sec) 
Latitude : 1959.40 minutes 
Longitude : 7034.40 minutes 

Ship GPS Antenna Height : -19.44 m 
Air Temperature : 15.30 C 

Relative Humidity : 69 % 

Air Pressure : 1016.40 mb 

Wind Speed : 4.0 m/s 

Sea Surface Temperature : 16.60 C 
Range : 373.21 m 

Theta : 1.71 rad 

Phi : 3.02 rad 

Heading : 68 degrees 

Bearing : 186.9 degrees 
Atmospheric Transmittance : 0.92 


Model's reference point is shown in the picture. This point is used to anchor the model to 


the ship image. It 1s located in x coordinate of the ship image area centroid and at the average y 
coordinate of the ship image's bottom line. 
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12. SHIP MODEL'S MAIN PLANES 


Description : 


This Appendix presents the position of the ship model's main planes. 


Conditions : 


Model elevation angles are 1.5 degree. The azimuth angles are 20 and 340 for the upper 
and lower pictures respectively. 


Notes : 


There are 4 main planes considered in the ship model. They are numbered from 0 to 3 and 
comprise the sides of the hull of the ship model. 
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13; MASK APPLICATION EXAMPLE 


Description : 


This Appendix presents an example ofthe mask application process corresponding to the 
plane 0 in the analysis box of the frame +0 of the unpolarized case of the BASESO.PTE file.There 
are 2 pictures in this Appendix. 


Conditions : 


Model elevation and azimuth angle are 2 and 70 degrees respectively. 
For the Image, from tables in Appendix E.2: 


File #50: 


Day  : April 9 1996 

Time : 62,187.5 seconds (elapsed time from 00h:00min local time — 17h31min8sec) 
Latitude : 1959.40 minutes 
Longitude : 7034.40 minutes 

Ship GPS Antenna Height : -19.44 m 
Aır Temperature : 15.30 C 

Relative Humidity : 69 % 

Aır Pressure : 1016.40 mb 

Wind Speed : 4.0 m/s 

Sea Surface Temperature : 16.60 C 
Range : 373.21 m 

Theta : 1.71 rad 

Phi : 3.02 rad 

Heading : 68 degrees 

Bearing : 186.9 degrees 
Atmospheric Transmittance : 0.92 


Notes : 


The upper picture is plotted just to show the position of the plane O in the image. The 
lower picture shows the real effect of the mask application in extracting the plane 0 corresponding 
pixels from the ship image. The small square in the picture is the reference point. It is shown in 
the pictures for better localization. The irregular left edge observed in the lower picture is the 
consequence of the ship mask application over the plane 0 mask. This application 1s necessary to 
avoid computation of the background part in the model plane 0 as if it was part of the ship for 
contrast calculation. Therefore, only the ship pixels inside the superimposed model's plane O are 
considered for contrast calculation. 
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APPENDIX I 
RESULT TABLES 


1. RESULT MATRIX DESCRIPTION 


Description : 


Table describing the columns of the result matrices. 


Conditions : 


None. 


Notes : 


None 
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RESULT MATRIX 


COLUMN DESCRIPTION 
DT : Elapsed Time 


Hotspot : Thermal Value 


~ 


о) 


Hotspot : x Position (data coordinates) 


Hotspot : y Position (data coordinates) 

Analysıs “Box” : Down Left Comer x Position (data coordinates) 
Analysıs “Box” : Down Left Comer y Position (data coordinates) 
Analysis “Box” : Up Right Corner x Position (data coordinates) 
Analysis “Box” : Up Right Comer y Position (data coordinates) 
Horizon Equation : Linear Coefficient (data coordinates) 

Horizon Equation : Angular Coefficient (data coordinates) 
Empty 

Empty 

Empty 

Empty 

Empty 

Empty 

Empty 

Empty 

Empty 


Empty 
Ship G.P.S. Position : Latitude (minutes) 


км 


к- 
= 


== 


= 
~) 


~ 


0 


~ 


] Ship G.P.S. Position : Longitude (minutes) 


һә 


2 Ship G.P.S. Position : Height (m) 


Aır Temperature (C) 


~ 
г. 


Relative Humidity (96) 
2 


in 


Aır Pressure (mb) 


bo 
ON 


Wind Speed (m/s) 


кә 
~J 


Sea Surface Temperature(C) 


t2 
00 


Camera Position : Latitude (minutes) 


== 
UJ 
© 


COLUMN DESCRIPTION 
Camera Position : Longitude (minutes) 


| 
3 
| 40 


Average Ship Stack Temperature (C) | 


29 
30 
] 


41 Ship Skin Temperature : point 1 (C) 


Ship Skin Temperature : point 2 (C) 


43 Ship Skin Temperature : point 3 (C) 


Ship Skin Temperature : point 4. (C) 


Ship Skin Temperature : point 5 (C) 







Ship Skin Temperature : point 6. (C) 
Ship Skin Temperature : point 7 (C) 










48 Ship Skin Temperature : point 8 (C) 
Ship Skin Temperature : point 9 (C) 





50 Ship Skin Temperature : point 10. (C) 


| Ship Skin Temperature : point 10 © _ 
















Ship Skin Temperature : point 15 (C) 
Ship Skin Temperature : point 16 (C) 
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COLUMN DESCRIPTION 
Contrast Ship x Background 


6 


pod 


Ship : Average Thermal Level 


ON 


2 Ship : Number of Pıxels 
Background : Average Thermal Level 


Background : Number of Pixels 


ON 
Un 


(ә 


Contrast Ship Above (above horizon line) x Sky 
Ship Above : Average Thermal Level 


ON 
~J 


Ship Above : Number of Pixels 


ON 
> 


Sky : Average Thermal Level 
Sky : Number of Pıxels 


~ 
e 


Contrast Ship Below (below horizon line) x Sea 
Ship Below : Average Thermal Level 


72 Ship Below : Number of Pixels 


~ 
UJ 


Sea : Average Thermal Level 


э 


Sea : Number of Pıxels 


“~ | ~ 
QU | 


Contrast Plane 0 x Background 


~J 


6 Plane 0 : Average Thermal level 


Plane 0 : Number of Pixels 


м 


7 


78 Plane 0 : Aspect Angle (degrees) 


э 
NO 


Contrast PlanelO x Background 
Plane 1 : Average Thermal level 


Plane ] : Number of Pixels 


00 


1 


00 


2 Рапе | : Aspect Angle (degrees) 


00 
95) 


Contrast Plane 2 x Background 
84 Plane 2 : Average Thermal level 


Plane 2 : Number of Pixels 


00 


5 


Plane 2 : Aspect Angle (degrees) 


00 
э 


Contrast Plane 3 x Background 


00 


8 Plane 3 : Average Thermal level 


Plane 3 : Number of Pixels 


[37 





%0 | Pime3 Aspect Angle (aegre _______________ 
Eco ы. 4. 
ж 
—> 













NOTES: 


1 . There is one Result Matrix for each polarization case. Therefore, since there are 3 
polarızatıon cases (Horizontal, Vertical and Unpolarized) in each BASEN.PTE file, there are 3 Result Matrices in each 
BASEN.PTE file as well. 

2 . The number of lines in each Result Matrix equals the number of frames of the 
correspondent polarization case. Therefore, each line corresponds to one frame data and thus, each Result Matrix is 100 
columns by “number of frames” lines . 

3 . Camera Referential is defined as follow: 


Ongin : Camera 
xaxis : True North 
yaxıs : True West 
zaxis : Upward 
4 . Contrast Values equal “ -1" means null average thermal level of the correspondent 
“target” part. Contrast Values equal “77" means null average thermal level of the correspondent “background” part. 
5 . Plane Aspect Angle is the angle between the plane normal vector and the line connecting 
the plane center and the camera. 
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Ze SUMMARY MATRIX DESCRIPTION 


Description : 
Table describing the columns of the “summary matrix”. 
Conditions : 
None. 


Notes: 


None 
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SUMMARY MATRIX 


COLUMN DESCRIPTION 


File Number (“N” in BASE”N”.PTE file) 
Average Ship Aspect Angle (degrees) 
Average Plane 0 Aspect Angle (degrees) 
Average Plane 1 Aspect Angle (degrees) 
Average Plane 2 Aspect Angle (degrees) 


Average Plane 3 Aspect Angle (degrees) 


Average Distance Ship - Camera (m) 


Ship Contrast Improvement : Horizontal Polarization 
Ship Contrast Improvement : Vertical Polarization 
Plane 0 Contrast Improvement : Horizontal Polarization 
Plane 0 Contrast Improvement : Vertical Polarization 


1 


— 


Plane 1 Contrast Improvement : Honzontal Polarization 


12 Plane 1 Contrast Improvement : Vertical Polarization 


= 


3 Plane 2 Contrast Improvement : Horizontal Polarızation 


Plane 2 Contrast Improvement : Vertical Polarization 


і 


5 Plane 3 Contrast Improvement : Horizontal Polarization 


р 


6 Plane 3 Contrast Improvement : Vertical Polarization 


Empty 
Empty 
Empty 


Hor. Polarization : Number of Frames 


= 
со 


2 


кә 
= 


Hor. Polarization : Average Ship Contrast 


~ 
м 


=3 


Hor. Polarization : Linear Coefficient for Ship Contrast Variation in Time 


N 
о 


Hor. Polarızation : Angular Coefficient for Ship Contrast Variation in Time 


bo 


00 


Hor. Polarization : Linear Correlation for Ship Contrast Variation in Time 


un 


2 Hor. Polarization : Standard Deviation for Ship Contrast Variation in Time 


2 Hor. Polarization : Average Plane 0 Contrast 
2 Hor. Polarization : Linear Coefficient for Plane O Contrast Variation in Time 


2 Hor. Polarization : Angular Coefficient for Plane 0 Contrast Variation in Time 
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COLU DESCRIPTION 


2 


чө 
NO 


Hor. Polarization : Linear Correlation for Plane O Contrast Variation in Time 


~ 


Hor. Polarization : Standard Deviation for Plane O Contrast Variation in Time 


UJ 
© 


1 Hor. Polarization : Average Plane ] Contrast 


32 Hor. Polarization : Linear Coefficient for Plane | Contrast Variation in Time 


UJ 


3 Hor. Polarization : Angular Coefficient for Plane 1 Contrast Variation in Time 


Hor. Polarization : Linear Correlation for Plane 1 Contrast Variation in Time 


UJ 
+ 


3 Hor. Polarization : Standard Deviation for Plane | Contrast Variation in Time 


5 
5 
40 


Hor. Polarızation : Standard Deviation for Plane 2 Contrast Variation in Time 


+ 
— 


Hor. Polarization : Average Plane 3 Contrast 


Hor. Polarization : Linear Coefficient for Plane 3 Contrast Variation in Time 


I 
кә 


A 


3 Hor. Polarization : Angular Coefficient for Plane 3 Contrast Variation in Time 


44 Hor. Polarization : Linear Correlation for Plane 3 Contrast Variation in Time 
45 Hor. Polarization : Standard Deviation for Plane 3 Contrast Variation in Time 
4 Ver. Polarization : Number of Frames 


Ver. Polarization : Average Ship Contrast 


00 


Ver. Polarızation : Linear Coefficient for Ship Contrast Variation in Time 


+ 
NO 


о 


Ver. Polanzation : Angular Coefficient for Ship Contrast Variation 1n Time 


Ver. Polarization : Linear Correlation for Ship Contrast Variation in Time 


Кл 
е 


Ver. Polarization : Standard Deviation for Ship Contrast Variation in Time 


in 
to 


Ver. Polarization : Average Plane 0 Contrast 


э 


Ver. Polarization : Linear Coefficient for Plane O Contrast Variation in Time 


un Ian ftn [tn [tn Un 
ON [| tA о 


4 


Ver. Polarization : Angular Coefficient for Plane O Contrast Variation in Time 
Ver. Polarization : Linear Correlation for Plane O Contrast Variation in Time 


Ver. Polarization : Standard Deviation for Plane O Contrast Variation in Time 


з 


Ver. Polarization : Average Plane | Contrast 


Ver. Polarization : Linear Coefficient for Plane | Contrast Variation in Time 


in 
00 
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COLUMN DESCRIPTION 
59 Ver. Polarization : Angular Coefficient for Plane 1 Contrast Variation in Time 
Ver. Polarization : Linear Correlation for Plane 1 Contrast Variation in Time 


Ver. Polarization : Standard Deviation for Plane 1 Contrast Variation in Time 


Ver. Polarization : Average Plane 2 Contrast 








B Ve ME 
вв | Ver Polarization i Linear Coeficient for Plane 3 Contrast Variation in Time | 
és | VenPotuizaion: Angular Coccen fr Plane 3 Contrast Variation in Time || 
70 | Ver Polaizain: Linear Correlation for Plane 3 Contrast Variation in Time || 
[п | Ver Polrizaon: Standard Deviation fr Plan 3 Contrast Variation in Time | - 
ee 


Unpolarized : Average Ship Contrast 


Е Unpolarized : Linear Coefficient for Ship Contrast Variation in Time 


Unpolarized : Angular Coefficient for Ship Contrast Variation in Time 





; 
Sa 










Unpolarized : Linear Correlation for Ship Contrast Variation in Time 


"XB Unpolarızed : Standard Deviation for Ship Contrast Variation 1n Time 
Unpolarized : Average Plane 0 Contrast 


Unpolarized : Linear Coefficient for Plane 0 Contrast Variation in Time 


тиге кеки ae 
0 i ne i i ot Vang OOO 
®© __| Unpotaried Standard Deviation for Plane O Contras Variation in Time 
IN со ағ 
ss | Unpolarized : Linear Correlation for Plane | Contrast Variation in Time | 
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COLUMN DESCRIPTION 
Unpolarized : Linear Coefficient for Plane 2 Contrast Variation in Time 
Unpolarızed : Angular Coefficient for Plane 2 Contrast Variation in Time 


Unpolarized : Linear Correlation for Plane 2 Contrast Variation in Time 


Unpolarized : Standard Deviation for Plane 2 Contrast Variation in Time 


ХО 


3 Unpolarızed : Average Plane 3 Contrast 


ХО 
> 


bo 


Unpolarized : Linear Coefficient for Plane 3 Contrast Variation in Time 


NO 


5 Unpolarized : Angular Coefficient for Plane 3 Contrast Variation in Time 


Unpolarized : Linear Correlation for Plane 3 Contrast Variation in Time 


хо 
~ 


Unpolarized : Standard Deviation for Plane 3 Contrast Variation in Time 
Ship : Degree of Polarization 

Plane 0: Degree of Polarization 

Plane 1: Degree of Polarization 

Plane 2: Degree of Polarization 


102 Plane 3: Degree of Polarization 


103 Ship / Honzontal Polarization : Average Thermal Level 
104 Ship / Vertical Polarization : Average Thermal Level 


105 Ship / Unpolarized : Average Thermal Level 


Plane 0 / Horizontal Polarization : Average Thermal Level 
107 Plane 0 / Vertical Polarization : Average Thermal Level 
Plane 0 / Unpolanzed : Average Thermal Level 

Plane 1 / Horizontal Polarization : Average Thermal Level 
Plane 1 / Vertical Polarization : Average Thermal Level 


Plane 1 / Unpolanzed : Average Thermal Level 


112 Plane 2 / Horizontal Polarization : Average Thermal Level 
113 Plane 2 / Vertical Polarization : Average Thermal Level 
114 Plane 2 / Unpolarızed : Average Thermal Level 

15 Plane 3 / Horizontal Polarızation : Average Thermal Level 


Plane 3 / Vertical Polarization : Average Thermal Level 


re 


li Plane 3 / Unpolarızed : Average Thermal Level 


Sea : Degree of Polarization 
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NOTES: 


1. Dimensions of Summary Matrix: 
130 columns(summary parameters) by 70 lines (BASEN.PTE files) 
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3. CONTRAST VARIATION IN TIME 


Description : 

This Appendix presents 3 tables with numerical results from the curve fitting process for 
the ship contrast variation in time for each of the 70 analyzed files. Each table corresponds to one 
polarization case (horizontal, vertical and unpolarized). The tables’ columns represent the 
calculated parameters and their lines represent the data files. 

Conditions : 

Appendix E.2 presents the conditions for each one of the 70 analyzed data files. 

Notes : 


None 
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SHIP CONTRAST VARIATION IN TIME 
CURVE FITTING DATA 
HORIZONTAL POLARIZATION CASE 


Contrast Coefficient Coefficient Coefficient Deviation 


6 |1246 (|14355 1-0.0418305 (|09803 101949072 | 
в 1163538 [161315 _ [ 1000563710 (|05569 = 1018452 | 
9 |21745 |20951 10.0194168 —  |0.762126 10.178305 —— 











139 |21021 |21562 1-0.0141800 10.600514 |05123 | 
-0.00548300 
37 |2.23930 
38 3.34082 
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0.386256 
0.225648 
1.30717 

0.654553 


File Average Linear Angular Correlation Standard 
Contrast Coefficient Coefficient Coefficient Deviation 


60 (|42967 |4161 10.0337203 [0.805949 |05626 - 
[66 3.52566 |35232 | -0.000940500 | 0.503029 10.741522 — | 
[68 [415409 4.04489 0.0258802_ 10.775146 10.241804 | 
[69 13.28224 13.0392  |005461 [0.790860 |05052 - 
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SHIP CONTRAST VARIATION IN TIME 
CURVE FITTING DATA 
VERTICAL POLARIZATION CASE 


Contrast Coefficient Coefficient Coefficient Deviation 
| 
5 |11404. |11979 -0.0118561 00 10.891219  |0.0835699 | 
6 [0.875345 ___| 0.966203 -0.0152129 |0.927931 |0.139760 — | 
|8 [1.45242  |154902 -0.0244514  — |0.822164 [0.183399 | 
9 . |17049. |12401 7 |012553 7 0942385  |0.648509 ___| 
Юи 10.754552, _ _ [0.716056 _ 5.|10:00608100. 10.579543 0.294911 ~ | 
0.169178 | 
[29 1212061 |1.96205 7 |0.0295708 [0.944763 [0.220237 | 
|36 |22480 |19706 7 [|0.0739158 | — (|0.998733 [0.340360 — | 

|37 |18554 [1.94024 |-0.0213637 [0.722923 10.228740 —— 
1.53497 
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сл 








Linear Correlation Standard 
Contrast Coefficient Coefficient Coefficient Deviation 
0.296846 
0.257858 
0.431382 
0.316965 
0.334339 
0.351712 
0.454030 
0.225624 
60 — 2.00306 |2.00704 -0.00120140 . [0.507623 [0.321723 | 
|66 [2.3059 |17801 ^ — 0.145194 — — |0.7726091] |11781, 
|68 3.80472 — |37336 7 |02018080 7 |0.65979 0.271108 - 
|69 |1.82662 7 |1.2004 -0.0222314  |077747  |0189958 | 
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SHIP CONTRAST VARIATION IN TIME 
CURVE FITTING DATA 
UNPOLARIZED CASE 


Contrast Coefficient Coefficient Coefficient Deviation 
6 [2.17321 |21458 7 |0.0064590 10591816 [0180736 | 


[8 |22559 |22560 0.000764200_ —  |0.507833 0.234531 | 
9. 122140 ]|2.32928 |-0.0147498 |09870 |00607 __ 


16 |10840 |1221 -0.0372685 |09570  |015964 | 
[18 — 0.798607 — |07081  |0023320 |07467 |020517 | 
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Contrast Coefficient Coefficient Coefficient Deviation 


[G0 |41455 |43462  |-0050645 |08346 [0326072 | 
|66 14.668556 |47919 |-0.0314910 |07875  |0.26526 | 
0.241617 
68 [551555 |57253 -0.0565878 [0.838031  0.585874 
C 129460 12759  |005605-7 |07148 | 0.386373 
0.266740 
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4. CONTRAST IMPROVEMENT AND DEGREE OF POLARIZATION 


Description : 

This Appendix presents a table with numerical results for ship to background contrast 
improvement and degree of polarization for each of the 70 analyzed files. The tables’ columns 
represent the calculated parameters and their lines represent the corresponding data files. 
Conditions : 

Appendix E.2 presents the conditions for each one of the 70 analyzed data files. 

Notes : 


None 
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SHIP CONTRAST IMPROVEMENT 


AND 
DEGREE OF POLARIZATION 
File Range Aspect Angle Horizontal Vertical Degree of 
(m) degree Improvement Improvement 
СЛ |998564 |57603 |1456 11 ] 0.850403 — à 1]-0.0184325 — | 


[3 199875 |16231 |09099 [0.618558  |-0018540 -Ҙ 
14 1186859 122280 | 0.898121 [0.713293  |-000045900 | 
[6 |182153 |72000  |05%6072 |04079  |00205949 | 
8 12413555 |99600: 107400 |06407 -0.00829560_ - 
[9 |241358 |33280 |09738 |07550 -0.0255484 | 
10 |24158 |54341 0.715161 106250 |00201 | 


RIO 


у 


[16 1336922 |72148 |1264 |10002 |-0.00250100 &—— 
18 15138.79 |13246 1089770 |1041  |0.0364165 | 
19 |50712 |82264 |17448 |10055 -0.0139722 __ 
0.834405 
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File Range Aspect Angle Horizontal Vertical Degree of 
Contrast Contrast Polarization 
m) (degree) Improvement Improvement 


1.14306 
0.879074 









[60 — | 323.227 | 361.143 11059 —1 à [0.483183 — [000741130 | 
66 |45945: |1570 |07512 10.493930 — |00850 - 
[68 1459698 | 161.190 |07510  |068%Я7 |-0.0214548 ____ 
69 |13950 |16287 11140 |06001  |-00702 - 





998.753 362.122 1.39508 0.608178 0.00218570 
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APPENDIX J 
RESULT GRAPHICS 


1. CONTRAST VARIATION IN TIME : DATA FILE SAMPLE 


Description : 


This Appendix presents a sample of the curve fitting process applied to the ship contrast 
variation in time. The 3 polarization cases are represented by corresponding plots. A comparative 
plot (with the 3 cases) is also presented. The data file used to generate this example was the 
BASEO1.PTE (file #01). 


Conditions : 
From the first line of the table in Appendix E.2: 
File 201 


Day : April 9 

Time: 44,018.0 seconds (elapsed time from 00h:00:min of April 9) 
Latitude : 1959.20 minutes 
Longitude : 7034.80 minutes 

Ship GPS Antenna height : -27.01 m 
Air Temperature : 15.30 C 

Relative Humidity : 72 % 

Aır Pressure : 1018.25 mb 

Wind Speed : 6.5 m/s 

Sea Surface Temperature : 18.55 C 
Range : 998.52 m 

Theta : 1.62 rad 

Fi : 2.41 rad 

Atmospheric Transmittance : 0.82 


Notes : 


None 
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SHIP CONTRAST VARIATION IN TIME : FILE SAMPLE ( BASEO1.PTE) 


HORIZONTAL POLARIZATION VERTICAL POLARIZATION 


3.5- Line Fitting : C = -0.021268 t + 2.1454 





3 
Uv 17 
2 (o 
E E 2.5 
О О Er ie, ри 
Line Fitting : C = -0.01597 t + 3.1053 1.5 
Mean = 2.0653 
Mean = 3.0426 1 
Std 2 0.11718 
Std = 0.2267 0.5 
0 
mr E). сы? 0 2 4 6 
Time (S) | Time (s) 
UNPOLARIZED CONTRAST : TIME VARIATION 
5 
4.5 
4 
3:5 
3 A 
% E қ | | | 
Б E25 | | | | 
2 УУ“ An МА 
Mean = 2.445 15 Horizontal Polarization 
Std 2 0.49104 1r ... Vertical Polarization 


0.5- —-Unpolarized 





0 2? 4 6 
Time (s) Time(s) 


2. CONTRAST VARIATION IN TIME : STATISTICAL DISTRIBUTION 


Description : 
This Appendix presents the statistical distribution (histograms) for the curve fitting 
parameters of the ship contrast variation in time. The 3 polarization cases are presented 


separately. All 70 data files were used to generate the plots in this Appendix. 


Conditions : 


Appendix E.2 presents the conditions for each one of the 70 files used in this Appendix. 


Notes : 
It can be observed by these histograms that the linear fitting parameters have a narrow 
distribution for all polarization cases. Therefore, they can be very well represented by their mean 


values. 
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3. CONTRAST IMPROVEMENT AND DEGREE OF POLARIZATION 


Description : 


This Appendix presents a the statistical distribution (histograms) for the ship contrast 
improvement and the image elements (ship, sea, sky) degree of polarization. The ship contrast 
improvement for horizontal and vertical polarization cases are presented in the same scale to 
permit easy comparison. The degree of polarization histograms for the image elements are also 
presented using the same scale to permit easy comparison. All 70 files data were used to generate 
the plots in this Appendix. 


Conditions : 


Appendix E.2 presents the conditions for each one of the 70 files used in this Appendix. 


Notes : 


The ship contrast improvement histograms show that in average, there is an advantage in 
using the horizontal filter for the scenario analyzed. On the other hand, the vertical filter results an 
average loss in contrast. The ship degree of polarization has a very narrow distribution when 
compared with the sea and the sky distributions. The sea exhibits predominance of vertical 
polarization. The ship degree of polarization is effectively zero and the sky shows a random 
behavior. 


157 


Егедиепсу 


Frequency 





SHIP CONTRAST IMPROVEMENT 


HORIZONTAL POLARIZATION 





Contrast Improvement 


VERTICAL POLARIZATION 





0.6 0.8 1 12 
Contrast Improvement _ . . /— 1 


158 


E. A чаиуиеутьтнти 











‚81126 
4209 


AN 


------- 





SHIP DEGREE OF POLARIZATION 


1 | 
Mean = -0.003478 
Std = 0.029823 


, 


Frequency 





Degree of Polarization 


SEA DEGREE OF POLARIZATION 





Degree of Polarization 


SKY DEGREE OF POLARIZATION 


Mean 7 -0.029317 
Std = 0.2416 
1 -0. 


8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1 
Degree of Polarization 







Frequency 






-------+-------:---+ nn s 


192 








4. GENERATION OF UNPOLARIZED IMAGES 


Description : 

This Appendix presents the analysıs plots for verifying the assumption of generating 
unpolarized images by averaging horizontally and vertically polarized images. There are three 
figures corresponding to the image's main elements: ship, sea and sky. In each picture, there are 
two plots: line fitting process and difference histogram. All 70 files data were used to generate 


the plots in this Appendix. 


Conditions : 


Appendix E.2 presents the conditions for each one of the 70 files used in this Appendix. 


Notes : 


None 
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SEA THERMAL LEVEL : COMPARISON UNPOLARIZED x AVERAGED “(V+H)/2” 
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SKY THERMAL LEVEL : COMPARISON UNPOLARIZED x AVERAGED “(V+H)/2” 
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S: TARGET ASPECT ANGLE EFFECTS 


Description : 

This Appendix presents the analysis plots for verifying the effect of the target aspect angle 
on its degree of polarization. There are six figures showing the degree of polarization of the ship 
and its main planes against the aspect angle at different ranges (distance sensor to target) For each 
plot, the data range information is indicated in its title. When “all ranges” is indicated, that means 
that all 70 files data were used to generate the plot. For the other plots (where a range interval is 


indicated) only the corresponding files were used. 


Conditions : 


Appendix E.2 presents the conditions for each one of the 70 files used in this Appendix. 


Notes : 
Column “range” in the table presented in Appendix E.2 corresponds to the range 


information used in the plots presented here (in this Appendix) 
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